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Abstract 


The  relationships  between  various  normalized  higher-order  concentration  moments  have 
been  investigated  using  a  large  data  set  of  concentration  fluctuations  obtained  in  a  boundary- 
layer  water  channel  with  high-resolution  laser-induced  fluorescence.  This  data  set  corre¬ 
sponds  to  a  series  of  comprehensive  measurements  of  plume  dispersion  in  a  number  of 
obstacle  arrays  (e.g.,  various  arrays  of  cubical  and  non-cubical  obstacles  in  aligned  and 
staggered  arrangements  with  uniform  and  random  heights).  A  remarkably  robust  feature 
of  all  the  concentration  data  was  the  observed  collapse,  of  the  third-  and  fourth-order 
normalized  concentration  moments  on  the  second-order  normalized  concentration  moment 
and  of  the  concentration  kurtosis  on  the  concentration  skewness,  to  a  series  of  “universal” 
curves.  These  “universal”  curves  were  identical  to  those  observed  previously  for  open-terrain 
plumes,  and  are  well  modeled  using  either  a  clipped-gamma  probability  density  function 
(PDF)  or  the  simpler  intermittent  exponential  PDF  for  the  concentration.  A  comparison 
of  the  shape  of  the  model  probability  distributions  to  the  measured  concentration  data  at 
various  plume  locations  showed  that  the  clipped-gamma  distribution  provided  a  good  rep¬ 
resentation  for  the  general  distribution  shape,  whereas  the  simpler  intermittent  exponential 
distribution  yielded  a  poor  conformance  to  the  measured  concentration  probability  distri¬ 
bution  (in  spite  of  the  fact  that  both  of  the  model  distributions  gave  a  good  representation 
for  the  first  four  concentration  moments).  The  implication  of  the  form  of  the  concentra¬ 
tion  PDF,  for  the  formulation  of  a  probabilistic  model  for  the  response  of  a  chemical  agent 
detector,  is  investigated  in  the  context  of  the  source  reconstruction  problem. 
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Resume 


On  a  etudie  les  relations  entre  des  moments  varies  de  concentrations  d’ordre  eleve  nor- 
malisees  en  utilisant  un  ensemble  important  de  donnees  de  fluctuations  de  concentrations 
obtenues  par  fluorescence  induite  au  laser  de  haute  resolution  dans  une  voie  d’eau  a  couche 
limite.  Get  ensemble  de  donnees  correspond  a  une  serie  de  mesures  comprehensives  de  dis¬ 
persions  en  panaches  dans  un  certain  nombre  de  reseaux  d’obstacles  (ex.  :  des  reseaux 
varies  d’obstacles  cubiques  et  non  cubiques,  arranges  en  lignes  ou  decales  et  de  hauteurs 
uniformes  ou  aleatoires).  On  a  observe  une  caracteristique  remarquablement  robuste  de 
toutes  les  donnees  de  concentrations,  consistant  en  I’affaissement  observe  des  moments  de 
concentrations  normalisees  de  troisieme  et  quatrieme  ordre  sur  le  moment  de  concentration 
normalisee  de  deuxieme  ordre  et  de  I’aplatissement  de  la  concentration  sur  I’etalement  de  la 
concentration,  en  une  serie  de  courbes  “  universelles  ” .  Ces  courbes  “  universelles  ”  etaient 
identiques  a  celles  observees  auparavant  pour  des  panaches  en  terrains  ouverts  et  elles  sont 
bien  modelisees  en  utilisant  soit  une  densite  de  probabilite  gamma  partielle  soit  une  densite 
de  probabilite  exponentielle  intermittente  plus  simple  de  la  concentration.  On  a  compare 
la  forme  des  distributions  de  probabilites  du  modele  avec  les  donnees  de  la  concentration 
mesuree  en  des  endroits  varies  du  panache  et  on  a  trouve  que  la  distribution  gamma  par¬ 
tielle  procurait  une  bonne  representation  de  la  forme  generale  de  la  distribution  alors  que  la 
distribution  exponentielle  intermittente  plus  simple  se  conformait  peu  a  la  distribution  de 
probabilites  de  la  concentration  mesuree  (en  depit  du  fait  que  les  deux  modeles  de  distribu¬ 
tion  donnaient  une  bonne  representation  des  quatre  premiers  moments  de  concentrations). 
On  etudie  actuellement  quelles  sont  les  implications  de  la  forme  de  la  densite  de  probabilite 
de  la  concentration  dans  la  formulation  du  modele  probabiliste  concernant  la  reponse  du 
detecteur  d’agents  chimiques,  dans  le  contexte  du  probleme  de  la  reconstruction  de  source. 
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Executive  summary 


The  Concentration  Probability  Density  Function  With  Implications 
for  Probabilistic  Modeling  of  Chemical  Warfare  Agent  Detector 
Responses  for  Source  Reconstruction 

E.  Yee;  DRDC  Suffield  TR  2008-077;  Defence  R&D  Canada  -  Suffield;  May  2008. 

Background:  A  critical  capability  gap  in  current  emergency  and  retrospective  manage¬ 
ment  efforts  directed  at  terrorist  incidents  involving  the  covert  release  of  chemical,  biologi¬ 
cal  and  radiological  (CBR)  agents  into  the  atmosphere,  is  the  determination  of  the  number 
of  sources  and  for  each  of  these  sources  the  location,  amount  of  agent  released,  and  time  of 
release  following  event  detection  by  a  network  of  CBR  sensors.  In  order  to  address  this  ca¬ 
pability  gap,  a  multi-national  collaborative  program  for  the  development  of  methodologies 
for  source  reconstruction,  involving  the  fusion  of  CBR  concentration  measurements  from 
remote  and  deployable  networks  of  sensors  with  model  concentration  data  obtained  from 
advanced  atmospheric  dispersion  models,  has  been  included  as  a  task  under  The  Technical 
Cooperation  Program  (TTCP),  CBR  Group  Technical  Panel  9  on  Hazard  Assessment. 

Principal  results:  A  probabilistic  model  for  the  response  of  a  chemical  agent  detector, 
that  provides  concentration  measurements  in  the  form  of  a  finite  number  of  discrete  bar 
levels,  is  formulated.  This  requires  that  the  probability  law  governing  the  instantaneous 
concentration  in  hazardous  plumes  dispersing  over  level,  unobstructed  terrain  and  built- 
up  (urban)  terrain  be  known.  Towards  this  purpose,  a  simple  probability  model  for  the 
plume  concentration  has  be  formulated  and  validated  using  a  large  data  set  of  measure¬ 
ments  of  plume  dispersion  in  idealized  obstacle  arrays.  This  data  set  was  obtained  from  a 
comprehensive  laboratary  study  in  a  water  channel  facility. 

Significance  of  results:  Many  operational  chemical  agent  detectors  provide  only  measure¬ 
ments  of  concentration  at  a  finite  number  of  discrete  levels  or  bars  (such  as,  for  example, 
the  Chemical  Agent  Monitor  (CAM)  developed  by  Smiths  Detection  or  the  AP4C  vapor  de¬ 
tector  developed  by  Proengin  Inc.).  The  probabilistic  model  for  a  bar- level  chemical  agent 
detector  developed  in  this  report  can  be  used  directly  in  source  reconstruction  method¬ 
ologies  that  have  been  developed  previously  to  determine  unknown  source  characteristics 
(e.g.,  location,  source  strength,  time  of  release,  etc.)  when  given  a  finite  number  of  noisy 
and  continuous  level  concentration  measurements  made  by  an  array  of  detectors.  With 
this  development,  the  sensor-driven  modeling  paradigm  can  be  integrated  into  operational 
warning  and  reporting  (information)  systems  that  combine  automated  data  acquisition, 
analysis,  source  reconstruction,  display,  and  distribution  of  CBR  hazard  prediction  and 
decision-support  products.  This  will  lead  to  a  more  complete  situational  awareness  in  cur¬ 
rent  warning  and  reporting  systems  for  the  operational  CBR  environment. 
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Future  work:  Future  plans  will  include  incorporating  the  proposed  probabilistic  model  for 
chemical  agent  bar-level  detectors  into  source  reconstruction  algorithms.  This  will  provide 
a  methodology  for  the  fusion  of  operational  chemical  agent  detector  data  (e.g.,  such  as  those 
obtained  from  in-service  CAM  or  AP4C  detectors)  with  model  predictions  of  chemical  agent 
dispersal  in  the  atmosphere. 
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Contexte  :  II  existe  une  lacune  importante  en  termes  de  capacitfe  en  efforts  de  gestion 
retroactive  des  situations  d’urgence  actuelles  comportant  I’mission  cachee  d’agents  chi- 
miques,  biologiques  et  radiologiques  (CBR)  dans  Tatmosphere  puisqu’il  faudrait  etre  capable 
de  determiner  le  nombre  de  sources,  et  pour  chacune  de  ces  sources,  le  lieu  d’emission,  la 
quantite  d’agent  emis  et  I’heure  de  remission  detectee  par  un  reseau  de  capteurs  CBR.  Pour 
etre  en  mesure  de  combler  cette  lacune,  le  Programme  de  Cooperation  technique.  Panel  9  du 
Croupe  technique  CBR  d’evaluation  des  dangers  comporte  une  nouvelle  tache.  Cette  tache 
consiste  en  un  programme  multinational  concerte  de  developpement  des  methodologies  de 
reconstruction  de  sources,  comprenant  notamment  la  fusion  des  mesures  des  concentra¬ 
tions  CBR  provenant  de  reseaux  de  capteurs  eloignes  pouvant  etre  deployes  et  ayant  des 
donnees  de  modeles  de  concentrations  obtenues  par  des  modeles  perfectionnes  de  dispersion 
atmospherique. 

Resultats  principaux  :  On  a  formule  un  modele  probabiliste  d’intervention  consistant 
en  un  detecteur  d’agent  chimique  qui  fournit  les  mesures  de  concentration  sous  forme  d’un 
nombre  fini  de  niveaux  de  barres  discrets.  Ceci  exige  de  connaitre  les  lois  de  probabilite  qui 
gouvernent  les  concentrations  instantanees  de  panaches  dangereux  se  dispersant  sur  un  ter¬ 
rain  nivele,  non  obstrue  et  construit  (urbain).  A  cette  fin,  on  a  formule  et  valide  un  modele 
simple  de  probabilites  de  la  concentration  des  panaches,  en  utilisant  un  ensemble  important 
de  donnees  des  mesures  de  la  dispersion  du  panache,  dans  des  echantillons  idealises  d’obs- 
tacles.  On  a  obtenu  cet  ensemble  de  donees  au  moyen  d’une  etude  comprehensive  dans  une 
installation  de  voie  d’eau  en  laboratoire. 

Portee  des  resultats  :  Beaucoup  de  detecteurs  d’agents  chimiques  operationnels  ne  four- 
nissent  que  les  mesures  de  concentration  a  un  nombre  fini  de  niveaux  ou  de  barres  discrets 
(tels  que  par  exemple  le  moniteur  d’agents  chimiques  -  CAM  -  mis  au  point  par  Smiths 
Detection  ou  le  detecteur  de  vapeurs  AP4C  mis  au  point  par  Proengin  Inc.).  Le  modele 
probabiliste  de  detecteur  d’agents  chimiques,  mis  au  point  dans  ce  rapport,  pent  etre  di- 
rectement  utilise  dans  les  methodologies  de  reconstruction  des  sources  developpees  aupa- 
ravant  pour  determiner  les  caracteristiques  des  sources  inconnues  (ex.  :  lieu,  puissance  des 
sources,  heure  d’emission,  etc.),  ceci  etant  donne  un  nombre  fini  des  mesures  des  niveaux 
de  concentration  bruyants  et  continus  effectuees  par  un  reseau  de  detecteurs.  Cette  mise  au 
point  permet  d’integrer  le  paradigme  de  modelisation  a  base  de  capteurs  dans  des  systemes 
operationnels  de  signal  d’avertissement  et  de  transmission  de  donnees  (informations)  qui 
combine  les  donnees  d’acquisition  informatisees,  les  analyses,  les  sources  de  reconstruction, 
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d’affichage  et  de  distribution  de  la  prediction  des  dangers  CBR  et  des  produits  de  soutien 
de  prise  de  decision.  Ceci  amenera  a  une  connaissance  plus  complete  de  la  situation  de 
I’environnement  operationnel  CBR  provenant  des  systemes  actuels  d’avertissements  et  de 
transmission  de  donnees. 

Perspectives  d’avenir  :  On  prevoit  d’inclure,  dans  les  algorithmes  de  reconstruction  de 
sources,  I’incorporation  de  ce  modele  probabiliste,  propose  ici,  aux  detecteurs  a  barres  de 
niveaux  d’agents  chimiques.  Ceci  devrait  fournir  une  methodologie  permettant  de  fusionner 
les  donnees  operationnelles  des  detecteurs  d’agents  chimiques  (comme  par  exemple  cedes 
obtenues  du  service  interne  de  detecteurs  CAM  ou  AP4C)  avec  le  modele  de  prediction  de 
la  dispersion  d’agents  chimiques  dans  I’atmosphere. 


VI 
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k,  s,  A  parameters  defining  clipped-gamma  distribution 

kv  von  Karman’s  constant  (~  0.4) 

C{bi\S)  likelihood  that  detector  reports  bar  level  bi  for  given  source  distribution  S 

Nf,  number  of  bar  levels  provided  by  a  chemical  agent  detector 

p{bi)  probability  for  detector  to  report  bar  level  bi 

Tc  integral  time  scale  for  concentration  fluctuations 

u*  friction  velocity  (=  (— u'rc')^/^) 

u'w'  Reynolds  shear  stress 

zq  roughness  length 

r(x)  gamma  function 

r(iy;x)  complementary  incomplete  gamma  function 

7  intermittency  factor  (probability  of  observing  a  non-zero  concentration  x) 

S(x)  Dirac  delta  function 

Af  frontal  area  index 

Ap  plan  area  index 

cji  measurement  noise  standard  deviation  for  concentration 

CJ2  standard  deviation  (uncertainty)  in  prediction  of  mean  concentration,  Cp 

CJ3  standard  deviation  (uncertainty)  in  prediction  of  mean-square  concentration,  c?p 

Ta  averaging  time  imposed  by  chemical  agent  detector 

X  instantaneous  (random)  concentration 

(•)  ensemble-averaging  operation 
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1  Introduction 


The  concentration  in  clouds  and  plumes  of  noxious  pollutants  (e.g.,  chemical  and  biologi¬ 
cal  warfare  agents)  is  inherently  a  random  variable,  owing  to  the  stochastic  nature  of  the 
turbulent  processes  responsible  for  dispersion  in  the  atmosphere.  The  often  neglected  prob¬ 
lem  of  concentration  fluctuations  in  turbulent  diffusion  has  received  increasing  interest  over 
the  past  25  years,  owing  to  the  fact  that  a  knowledge  of  the  characteristics  of  concentra¬ 
tion  fluctuations  is  of  considerable  importance  in  many  industrial  and  environmental  fluid 
mechanics  problems.  For  example,  the  statistical  properties  of  concentration  fluctuations 
in  a  dispersing  plume  are  important  to  the  assessment  of  risk  from  the  release  of  certain 
highly  toxic  materials  (e.g.,  industrial  chemicals,  chemical  warfare  agents)  in  which  there  is 
a  nonlinear  relationship  between  concentration  and  duration  of  exposure  for  a  given  level 
of  harmful  effect  [1],  [2].  Similarly,  short-term  concentration  fluctuations  are  very  relevant 
to  estimating  ignition  hazard  from  the  leakage  of  flammable  gases  (e.g.,  fuel-air  mixtures, 
liquified  natural  gas  spills)  in  which  it  is  necessary  to  determine  the  probability  that  the 
instantaneous  concentration  lies  between  the  lower  and  upper  flammability  limits  [3].  Still 
other  important  applications  include  the  prediction  of  the  probability  of  visibility  through 
obscurant  clouds,  the  characterization  of  the  perception  of  odours  required  to  evaluate  the 
nuisance  due  to  malodourous  substances,  and  the  understanding  of  fast  nonlinear  phys- 
iochemical  processes  required  for  the  design  of  efficient  mixing  and  combustion  processes 

[4]. 

The  statistical  description  of  the  natural  random  fluctuations  in  the  instantaneous  concen¬ 
tration  of  a  plume  dispersing  in  the  atmosphere  is  conveniently  embodied  in  the  probability 
density  function  (PDF)  of  concentration.  More  specifically,  in  order  to  investigate  and 
model  concentration  fluctuation  phenomenology,  it  is  useful  to  know  the  probability  den¬ 
sity  function  of  concentration  which  necessarily  embodies  all  the  higher-order  concentration 
moments  of  the  stochastic  process.  Chatwin  [5]  discusses  the  importance  of  the  concentra¬ 
tion  probability  density  function  with  reference  to  the  assessment  of  hazards  posed  by  the 
release  of  flammable  and  toxic  gases.  In  particular,  he  emphasizes  the  need  to  characterize 
precisely  the  upper  tail  of  the  probability  density  function  in  order  to  determine  exceedances 
of  critical  threshold  levels  by  peak  concentrations. 

The  shape  and  form  of  the  concentration  PDF  of  dispersing  clouds  and  plumes  of  contam¬ 
inant  in  the  atmosphere  have  been  studied  by  numerous  investigators.  In  consequence,  a 
number  of  models  have  been  proposed  for  the  concentration  PDF.  The  lognormal  distribu¬ 
tion  for  concentration  was  proposed  by  Csanady  [6]  and  used  by  Jones  [7].  The  exponential 
distribution  was  advocated  by  Barry  [8]  and  applied  by  Hanna  [9].  The  clipped-normal 
distribution  was  proposed  by  Lewellen  and  Sykes  [10]  through  application  of  the  princi¬ 
ple  of  maximum  entropy  and  adopted  by  Sawford  [11]  and  Mylne  and  Mason  [12]  as  the 
distribution  that  best  represents  the  concentration  fluctuation  data.  Yee  et  al  [13]  tested 
seven  candidate  probability  distributions  including  the  gamma,  Weibull,  and  conjugate 
beta  distributions  with  reference  to  concentration  fluctuation  measurements  obtained  in 
a  water  channel.  Finally,  Yee  and  Chan  [14]  proposed  a  clipped-gamma  distribution  for 
the  concentration  PDF  and  demonstrated  that  this  distribution  provided  excellent  agree- 
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ment  with  the  observed  higher-order  concentration  moment  relationships  obtained  from 
high-resolution  concentration  fluctuation  measurements  made  during  the  cooperative  Con¬ 
centration  Fluctuation  Experiments  (CONFLUX)  project  [15],  [16]. 

All  the  research  effort  cited  above  for  the  modeling  of  the  underlying  distributional  shape 
of  the  fluctuating  concentration  focussed  exclusively  on  plumes  dispersing  over  level,  unob¬ 
structed  terrain.  However,  it  needs  to  be  emphasized  that  as  the  fraction  of  the  World’s 
population  that  lives  in  cities  continues  to  grow,  it  is  becoming  increasingly  important  to 
address  the  urgent  problem  of  modeling  the  dispersion  of  toxic  materials  in  the  urban  en¬ 
vironment.  This  modeling  is  challenging  because  it  is  necessarily  characterized  by  the  wide 
spectrum  of  length  and  time  scales  associated  with  the  interaction  of  the  wind  flow  with  the 
complex  geometries  and  interfaces  (building  configurations,  network  of  streets)  that  com¬ 
prise  the  urban  fabric.  Presently,  there  is  a  paucity  of  information  concerning  the  statistical 
description  of  concentration  fluctuations  in  plumes  dispersing  in  urban  (built-up)  areas  and, 
more  particularly,  the  characterization  of  the  shape  of  the  one-point  concentration  PDF  for 
instantaneous  plumes  dispersing  within  urban  areas. 

Nevertheless,  a  small  number  of  studies  have  focussed  on  the  elucidation  of  the  statistical 
characteristics  of  concentration  fluctuations  in  plumes  dispersing  in  an  urban  environment. 
For  example,  Yee  and  Biltoft  [17]  extracted  a  number  of  concentration  statistics  (e.g.,  con¬ 
centration  variance,  concentration  PDF,  various  concentration  time  and  length  scales  of 
dominant  plume  motions)  in  a  plume  dispersing  through  a  large  array  of  building-like  ob¬ 
stacles.  Gailis  and  Hill  [18]  reported  a  wide  range  of  concentration  statistics  and  other 
quantitative  descriptors  of  plume  behaviour  for  tracer  dispersion  within  a  large  array  of 
obstacles  in  a  boundary-layer  wind-tunnel  simulation.  Yee  et  al  [19]  provided  detailed  com¬ 
parisons  of  concentration  statistics  in  a  plume  dispersing  through  a  large  array  of  building¬ 
like  obstacles  at  three  different  scales;  namely,  at  full-scale  in  a  field  experiment,  at  1:50 
scale  in  a  wind-tunnel  simulation,  and  at  1:205  scale  in  a  water-channel  simulation.  Finally, 
Klein  et  al  [20]  analyzed  and  compared  concentration  fluctuation  measurements  from  the 
Joint  Urban  2003  (JU2003)  full-scale  and  wind-tunnel  experiments. 

The  objective  of  this  paper  is  to  investigate  how  the  structure  of  the  statistical  proper¬ 
ties  of  a  plume  dispersing  through  regular  arrays  of  obstacles,  as  manifested  through  the 
higher-order  moments  of  concentration,  are  modified  in  comparison  to  those  observed  for 
a  plume  dispersing  over  a  level,  unobstructed  terrain.  The  measured  relationships  between 
the  various  concentration  moments  for  the  array  plumes  will  be  compared  with  predictions 
of  these  relationships  obtained  from  some  simple  models  for  the  probability  density  function 
of  concentration  fluctuations.  The  utility  of  these  simple  probability  distribution  models  to 
characterize  the  shape  of  the  measured  concentration  probability  distribution  for  an  array 
plume  is  assessed.  The  implication  of  the  form  of  the  concentration  PDF  for  the  formula¬ 
tion  of  a  probabilistic  model  for  the  response  of  a  chemical  agent  detector,  that  provides 
concentration  measurements  only  at  a  number  of  discrete  levels  or  bars,  is  investigated  in 
the  context  of  the  source  reconstruction  problem. 
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2  Models  for  concentration  PDF 


For  a  fixed  point  in  space  x  and  a  given  instant  in  time  t,  the  random  concentration  field 
x(x,  t)  >  0  of  a  hazardous  pollutant  dispersing  in  the  atmosphere,  can  be  characterized  by 
the  one-point  concentration  PDF  /(c;x,t)  defined  by 


/(c;x,t)(ic  =  Pr{c  <  x(x,t)  <c  +  dc}. 


(1) 


Here,  the  right-hand  side  denotes  the  probability  that  the  random  variable  x(x,  t)  falls 
in  the  interval  of  sample  space  values  between  c  and  c  +  dc  for  different  realizations  of 
the  turbulent  dispersion  of  the  hazardous  cloud  or  plume.  In  Yee  and  Chan  [14],  a  new 
concentration  PDF  model  corresponding  to  a  left-shifted  and  clipped-gamma  distribution, 
was  proposed.  In  this  model,  the  concentration  PDF  has  the  form 


/(c;  x,t) 


C  +  XY  ^  exp(-(c  +  A)/s)  , 

— j  — sm — 


(2) 


with  k  =  k{'x.,t)  >  0,  s  =  s(x, t)  >  0,  A  =  A(x, t)  >  0,  7  =  7(x, t)  G  [0,1],  r(x)  is  the 
gamma  function,  S{x)  is  the  Dirac  delta  function  and  the  range  for  c  is  0  <  c  <  00.  The 
total  PDF  of  concentration  in  Eq.  (2)  is  composed  of  a  mixed  fluid  part  (first  term  on 
right-hand  side  of  equation)  that  results  from  in-plume  mixing  of  eddies  that  contain  the 
scalar  contaminant,  and  an  unmixed  ambient  fluid  part  (second  term  on  right-hand  side 
of  equation)  that  is  produced  by  plume  meandering  producing  intermittent  periods  of  zero 
concentration  for  a  fraction  of  time  (1  —  7).  Alternatively,  7  =  Pr{x(x,  t)}  >  0  is  the 
intermittency  factor  that  determines  the  probability  of  observing  a  non-zero  concentration 
X  at  a  point  (x,  t)  in  space-time. 


The  concentration  PDF  in  Eq.  (2)  is  completely  determined  by  four  parameters:  namely,  7, 
k,  s  and  A.  However,  for  practical  applications,  we  require  a  simple  concentration  PDF  that 
predicts  an  arbitrary  moment  of  the  concentration  given  information  on  only  the  two  lowest- 
order  moments  of  concentration;  namely,  the  mean  concentration  (first-order  moment)  and 
the  mean-square  concentration  (second-order  moment).  In  other  words,  we  require  a  simple 
PDF  form  that  can  be  specified  with  no  more  than  two  parameters,  and  yet  is  capable  of 
fitting  observed  cloud  or  plume  concentration  data  over  a  wide  range  of  atmospheric  con¬ 
ditions  (various  diabatic  conditions)  and  terrain  types  (level,  unobstructed  terrain,  urban 
terrain,  etc.).  This  constraint  arises  from  the  fact  that  currently,  the  most  advanced  disper¬ 
sion  models  have  difficulty  in  predicting  even  the  two  lowest-order  concentration  moments 
with  sufficient  accuracy  [21]. 


To  provide  a  special  form  of  the  clipped-gamma  distribution  that  is  determined  by  two 
(rather  than  four)  parameters,  Yee  and  Chan  [14]  introduced  two  additional  constraints. 
Firstly,  it  was  assumed  that  7  can  be  uniquely  determined  as  the  area  remaining  under  the 
gamma  PDF  curve  for  c  >  0  after  a  left-shift  of  c  by  the  amount  A,  so 


7  =  7(fc,s,  A) 


F(fe;A/s) 


exp(— c/s) 
sT{k) 


dc 


m 


(3) 
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where  T{v]x)  denotes  the  complementary  incomplete  gamma  function.  With  this  assump¬ 
tion,  we  have  7  ^  1“  (i.e.,  non-intermittent  plume  concentration)  when  A  ^  O"*"  (i.e.,  no 
left-shift)^.  Secondly,  concentration  data  measured  in  the  CONFLUX  project  [15],  [16]  was 
used  to  provide  a  simple  relation  between  the  normalized  mean-square  concentration  and 
plume  intermittency: 

((X/Cf)  =  3/7.  (4) 

where  (•)  denotes  the  ensemble  average  (viz.,  average  over  a  number  of  independent  real¬ 
izations  of  a  process)  and  C  is  the  ensemble-averaged  concentration  (viz.,  C  =  (x)). 


Note  that  with  the  parameterization  for  7  given  by  Eq.  (4),  the  plume  concentration  is 
implicitly  assumed  to  be  non-intermittent  (viz.,  7=1)  for  ((x/C*)^)  G  [1)3]^  because  7 
cannot  exceed  unity  in  value.  For  a  non-intermittent  plume  concentration,  the  concentration 
PDF  of  Eq.  (2)  reduces  to  a  simple  gamma  PDF  as  7  =  1  and  A  =  0  (which  occurs  when 
((x/C*)^)  G  [1,3]);  viz.,  the  concentration  PDF  reduces  to 


/(c;x,t) 


exp(— c/s) 
sF(A:)  ^ 


(5) 


when  ((x/C)^)  ^  [1,  3].  More  succinctly,  this  special  form  of  the  clipped-gamma  distribution 
is  defined  with  7  determined  as 


so  that  when  ((x/C')^)  £  [1,3],  7  =  1  and  A  =  0;  and,  the  clipped-gamma  distribution 
simplifies  to  the  simple  gamma  distribution  in  this  case. 


The  n-th  concentration  moment  of  Eq.  (2)  is  given  by 


(x") 


c’^/(c;  X,  t)  dc 


j  +  k]  A/s) 

r(fc) 


n  G  N. 


(7) 


With  the  constraints  given  by  Eqs.  (3)  and  (6),  the  model  concentration  PDF  parameters 
k,  s  and  A  in  Eq.  (2)  can  be  obtained  by  solving  the  following  system  of  transcendental 
equations; 


1 

s 


7  + 


r(A;) 


k 

exp(-A/s); 


(8) 


(^(A/s)7  T  ( — A/s  k  l)/s^ 

(-A/s  -I-  /c)7  -I-  (A/s)^  exp(-A/s)/r(A:)]  ^ 


(9) 


^Note  that  a:  — >  a“  and  a:  — >  denotes  the  approach  to  x  to  the  limiting  value  a  from  the  left  (viz., 

from  values  of  x  smaller  than  a)  and  from  the  right  (viz.,  from  values  of  x  larger  than  a),  respectively. 

^The  concentration  variance  (x^^)  =  (x^)  —  >  0  (x^  =  X  “  U),  so  (x^)  >  or,  equivalently, 

((x/D)^)  >  1  is  the  lower  bound  for  the  normalized  mean-square  concentration. 
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and 


7  =  mm 


(10) 


1,3 


r(fc;A/s) 

m 


Thence,  for  a  specified  value  of  {(x/C)'^),  Eqs.  (9)  and  (10)  can  be  solved  for  k  and  (A/s),^ 
and  these  values  can  be  subsequently  substituted  into  Eq.  (8)  to  obtain  s,  after  which  the 
value  of  A  can  be  deduced.  The  solution  of  these  equations  for  k,  s  and  A  as  a  function  of 
((x/C*)^)  is  exhibited  in  Eigure  1. 


<  ( j(/cf  > 

(a)  Clipped-gamma  parameter  A 


(b)  Clipped-gamma  parameter  k 


(c)  Clipped-gamma  parameter  s 


Figure  1:  Dependence  of  the  parameters  (k,  s,  and  \)  of  a  clipped-gamma  distribution  on 
the  normalized  mean-square  concentration  ((x/C)^). 


The  skewness  S  and  kurtosis  K  are  defined  by 


_  fa”)  ((x  -  C)3> 


iix-c)^) 

/t4  /t4 


(11) 


where  =  (x^^)  =  {{x  ~  C*)^)  is  the  variance  of  concentration.  The  skewness  and  kurtosis 
can  be  readily  computed  from  the  relationships  between  the  central  concentration  moments 

®Note  that  the  factor  1/s  in  Eq.  (9)  can  be  expressed  as  a  function  of  only  k  and  (A/s)  using  Eq.  (8). 
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(moments  about  the  mean  concentration)  and  the  concentration  moments  about  the  origin 
[cf.  Eq.  (7)].  It  can  be  easily  verified  that  the  second-,  third-,  and  fourth-order  central 
concentration  moments  are  related  to  the  concentration  moments  about  the  origin  as  follows: 


and 


al  ^  =  {X^)  -  C\  (12) 

(x'')  =  (x')-3C(x2)  +  2C'3,  (13) 


(X'")  =  (x")  -  4C(x3)  +  6C2(x2)  - 


(14) 


In  view  of  the  complexity  of  the  clipped-gamma  distribution  of  Eq.  (2) ,  it  is  useful  to  consider 
a  simpler  model  for  the  concentration  PDF;  namely,  the  exponential  PDF  for  intermittent 
concentrations  proposed  by  Barry  [8].  This  concentration  PDF  is  a  two-parameter  PDF 
that  has  the  following  simple  form; 

2 

/(c;  ^  (“^)  +  “  7)<J(c).  (15) 

For  this  distribution,  the  intermittency  factor  is  related  to  the  normalized  mean-square 
concentration  as  follows: 

Furthermore,  the  relationship  between  the  normalized  higher-order  and  second-order  con¬ 
centration  moments  for  the  intermittent  exponential  PDF  is  given  by 


( 


r(n+l) 

2^-1  \VC'/  / 


n  G  N. 


(17) 


Strictly,  the  exponential  PDF  imposes  the  constraint  that  ((x/C*)^)  ^  2  [cf.  Eq.  (16)],  owing 
to  the  fact  that  7  <  1.  However,  similar  to  the  clipped-gamma  distribution,  we  implicitly 
assume  that  7  =  1  when  ((x/C)^)  ^  [1)2].  In  this  case,  the  intermittent  exponential  PDF 
given  by  Eq.  (15)  reduces  to  the  usual  exponential  PDF;  viz., 

/(c;x,t)  =  ^exp  .  (18) 

To  summarize,  7  for  the  intermittent  exponential  distribution  is  defined  as 

so  when  ((x/C*)^)  £  [1)2],  7=1  and  the  intermittent  exponential  distribution  simplifies  to 
the  usual  exponential  distribution  of  Eq.  (18). 

It  is  interesting  to  note  that  the  intermittent  exponential  distribution  given  by  Eq.  (15)  can 
be  considered  to  be  a  special  case  of  the  clipped-gamma  distribution  given  by  Eq.  (2).  In 
this  form  of  the  clipped-gamma  distribution,  Eq.  (19)  is  used  to  relate  the  intermittency 
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factor  to  the  normalized  mean-square  concentration,  rather  than  the  ansatz  of  Eq.  (6). 
From  Eq.  (3),  this  immediately  implies  that 


Solving  Eqs.  (20)  and  (9)  for  k  and  A/s,  it  can  be  shown  that  k  =  1  and  7  =  exp(— A/s), 
which  by  virtue  of  Eq.  (8)  implies  that  s  =  I/7.  The  latter  value  for  s  is  appropriate  for 
the  normalized  concentration  y/C  (or  equivalently,  for  a  system  of  units  where  (7=1).  For 
the  unnormalized  concentration  7,  s  =  C /^.  Using  these  values  for  k,  7  and  s  in  Eq.  (2), 
the  concentration  PDF  for  this  special  case  reduces  to  the  form: 

/(c;x,t)  =  -exp(-(c/(7  +  A)/s) -h  (1  -  7)(5(c) 

=  ^exp  +  (1 -7)5(c),  (21) 

which  is  exactly  the  intermittent  exponential  PDF  of  Eq.  (15). 

3  Experimental  data  for  model  evaluation 


Previously,  a  comprehensive  concentration  fluctuation  data  set  obtained  in  field  experi¬ 
ments  during  the  September  1991,  November  1992,  May  1993,  and  May  1994  phases  of 
the  CONFLUX  project  was  used  for  studying  the  relationships  between  the  higher-order 
concentration  moments  [14]  under  a  very  wide  range  of  atmospheric  conditions.  In  this 
study,  1,107  individual  CONFLUX  concentration  time  series  from  a  large  number  of  differ¬ 
ent  experiments  covering  a  very  wide  range  of  atmospheric  conditions  were  used.  All  these 
experiments  involved  the  continuous  and  controlled  release  of  propylene  (CsHg)  tracer  gas 
into  the  atmosphere.  The  measurements  were  made  at  downwind  fetches  x  of  between  12.5 
and  1,000  m  from  the  source,  under  moderately  convective  to  extremely  stable  atmospheric 
conditions  in  which  mechanical  turbulence  was  suppressed  by  the  stable  stratification.  The 
concentration  time  series  were  measured  over  a  wide  range  of  receptor  positions  in  both 
lateral  and  vertical  cross-sections  through  the  plume  —  lateral  plume  positions  varied  from 
the  mean-plume  centerline  at  y/uj^  =  0  to  the  extreme  plume  fringes  at  yjcFy  ~  ±4,  where 
y  is  the  crosswind  distance  from  the  mean-plume  centerline  and  ay  is  the  mean-plume  dis¬ 
persion;  and,  vertical  positions  ranged  from  0.5  to  16  m  height  above  ground.  The  source 
was  placed  at  heights  ranging  from  1  to  5  m. 

The  analysis  conducted  by  Yee  and  Chan  [14],  using  this  concentration  fluctuation  data  set, 
showed  that  despite  the  diversity  of  the  different  field  experiments  which  were  undertaken 
under  a  wide  range  of  atmospheric  conditions,  scatterplots  of  various  higher-order  normal¬ 
ized  concentration  moments  {{x/C)^)  {n  =  3,  4,  5,  6,  7,  and  8)  versus  the  second-order 
normalized  concentration  moment  ((x/C)^)  each  exhibited  a  remarkable  collapse  of  the 
data  onto  a  single  curve.  This  collapse  is  a  remarkably  robust  feature  of  the  concentration 
data  and  provides  compelling  evidence  that  the  concentration  PDF  can  be  represented  (to  a 
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very  good  approximation)  with  2  parameters;  namely,  one  for  location  (mean  concentration) 
and  one  for  scale  (root  mean-square  concentration,  or  equivalently,  concentration  standard 
deviation).  Furthermore,  Yee  and  Chan  [14]  demonstrated  that  the  clipped-gamma  PDF 
provided  good  predictions  of  the  observed  relationships  between  various  higher-order  nor¬ 
malized  concentration  moments  and  the  second-order  normalized  concentration  moment. 

This  collapse  of  concentration  data  was  demonstrated  using  measurements  of  plume  disper¬ 
sion  over  level,  unobstructed  terrain.  In  view  of  this,  the  question  that  we  pose  is  as  follows. 
Do  the  concentrations  in  a  plume  dispersing  through  and  over  an  urban  environment  con¬ 
sisting  of  a  large  collection  of  buildings  and  other  obstacles  (e.g.,  cars  lining  a  street,  treed 
areas  in  city  green  spaces)  exhibit  a  collapse  of  the  higher-order  concentration  moments 
on  the  second-order  concentration  moment?  If  so,  can  this  collapse  be  predicted  using 
the  clipped-gamma  distribution  as  in  the  case  of  plume  dispersion  over  level,  unobstructed 
terrain? 

To  address  these  questions,  we  will  use  a  large  concentration  data  set  obtained  from  mea¬ 
surements  of  the  near-held  dispersion  of  tracer  plumes  in  different  arrays  of  building-like 
obstacles.  This  comprehensive,  high-quality  data  set  of  dispersion  of  array  plumes  was 
obtained  from  water  channel  simulations  carried  out  at  Coanda  Research  &  Development 
Corporation’s  (Burnaby,  British  Columbia,  Canada)  boundary-layer  water  channel  [22]. 
The  water  channel  is  specially  designed  for  dispersion  modeling  and  has  a  working  section 
of  10  m  length,  1.5  m  width,  and  0.9  m  height.  The  upwind  part  of  the  working  section  of 
the  water  channel  was  used  to  generate  a  naturally  grown  neutrally-stratihed  rough- walled 
boundary  layer  of  about  0.3  m  thickness.  This  was  accomplished  using  a  “turbulence”  grid 
made  of  square  bars  19  mm  x  19  mm  placed  at  the  start  of  the  channel  inlet  to  generate 
turbulence  at  the  entry  of  the  working  section.  A  sawtooth  fence  with  a  base  width  equal 
to  that  of  the  channel  and  a  height  of  70  mm  was  placed  200  mm  downstream  of  the  square 
bar  array.  The  combination  of  sawtooth  fence  and  turbulence  grid  was  used  to  trigger  and 
produce  some  of  the  larger-scale  eddies  in  the  boundary  layer. 

Downstream  of  the  sawtooth  fence,  the  floor  of  the  water  channel  was  covered  with  a  black 
anodized  expanded  metal  mesh  of  height  4  mm  with  a  total  streamwise  length  of  6  m  to 
give  a  sufficiently  long  fetch  of  uniform  surface  roughness  for  the  upstream  flow  to  develop. 
The  expanded  metal  mesh  had  diamond-shaped  holes  that  were  approximately  11  mm  long 
in  the  streamwise  direction  and  25  mm  wide  in  the  spanwise  direction.  After  this  uniform 
upstream  fetch  of  expanded  metal  mesh,  the  flow  then  encountered  the  model  obstacle  array. 
Downstream  of  the  obstacle  array,  the  flow  encountered  a  section  covered  with  the  same 
black  anodized  expanded  metal  mesh  that  was  used  for  the  initial  upstream  boundary-layer 
development. 

Measurements  of  the  vertical  profiles  of  mean  velocity,  u,  streamwise  velocity  variance, 
and  shear  stress  r  =  u'w'  made  using  a  4-beam,  2-component  fibre-optic  laser 
Doppler  velocimeter  (LDV)  within  the  equilibrium  boundary  layer  upstream  of  the  obstacle 
array,  showed  that  these  profiles  were  homogeneous  to  within  ±5%  in  the  spanwise  direction 
across  the  channel.^  The  boundary-layer  thickness,  6,  taken  to  be  the  height  where  the  mean 

^Overbars  and  primes  denote  time  averages  and  departures  therefrom;  and,  (u,  v,  w)  are  the  instantaneons 
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wind  speed  reached  99%  of  the  free  stream  value,  was  found  to  be  275  mm.  At  this  point, 
the  mean  wind  speed  us  was  0.375  m  s“^.  The  friction  velocity  u*  =  (— determined 
from  measurements  of  the  shear  stress  in  the  constant  stress  layer  near  the  surface  was 
0.0255  m  s“^,  giving  u^^/us  =  0.068  for  the  water  channel  simulations.  A  least-squares  fit  of 
the  usual  log-law  profile  for  the  mean  wind  speed  in  a  regular  rough-wall  boundary  layer, 
u/u^  =  \og\{z  —  d) / zo\ ,  where  zq  is  the  roughness  length,  d  is  the  zero-plane  displacement, 
and  kv  ~  0.4  is  von  Karman’s  constant  gave  the  following  results:  zq  =  0.35  ±  0.05  mm, 
assuming  a  zero-plane  displacement  d  of  2.8  mm  (using  the  common  rule  of  thumb  that  d 
should  be  approximately  70%  of  the  height  of  the  roughness  elements). 

Six  different  types  of  obstacle  arrays  were  used  in  this  study.  The  first  obstacle  array  is 
shown  in  Figure  2.  This  obstacle  array  corresponds  to  a  1:205  scale  of  a  large  array  of 
building-like  obstacles  used  for  the  Mock  Urban  Setting  Test  (MUST)  field  experiment 
conducted  in  September  2001  at  U.S.  Army  Dugway  Proving  Ground  [17].  For  the  water 
channel  simulation  of  the  MUST  field  experiment,  each  obstacle  was  a  rectangular  paral- 
lelipiped  with  a  width  of  59.4  mm,  length  of  11.8  mm,  and  height  of  12.4  mm.  A  total  of  120 
obstacles  was  placed  in  an  aligned  configuration  consisting  of  12  rows  of  10  obstacles  with 
an  obstacle  spacing  (face-to-face)  of  38.5  mm  in  the  spanwise  (y)  direction  and  62.9  mm  in 
the  streamwise  (x)  direction  (cf.  Figure  2).  The  frontal  (A/)  and  plan  (Ap)  area  indices  of 
the  obstacle  array  were  0.10  and  0.096,  respectively.®  The  MUST  obstacles  for  the  water 
channel  simulation  were  machined  from  anodized  aluminum  and  installed  on  a  lexan  panel. 

The  remaining  five  obstacle  arrays  were  composed  of  cubes  and/or  parallelipipeds  (rect¬ 
angular  blocks),  constructed  from  Lego"’"'^  blocks.  These  obstacle  arrays  were  constructed 
from  a  square  ‘repeating’  unit  cell  (63.5  mm  x  63.5  mm)  shown  in  Figure  3.  The  unit  cell 
was  divided  into  4  quadrants  labelled  A,  B,  C,  and  D.  Each  obstacle  (cube  or  rectangular 
block)  of  the  array  had  a  square  cross  section  with  a  side  length  H  =  31.75  mm.  An  obstacle 
can  be  placed  on  any  one  of  these  4  quadrants  in  the  unit  cell.  The  obstacle  arrays  were 
constructed  from  a  16  x  16  unit  cell  array  layout  placed  on  a  Lego"*"^  baseplate. 

Schematic  drawings  of  the  obstacle  arrays  constructed  using  cubes  and/or  parallelipipeds 
are  shown  in  Figures  4  to  8.  Each  of  these  arrays  had  a  plan  area  index  Ap  =  0.25.  Eor  the 
arrays  of  aligned  obstacles,  the  obstacles  were  placed  in  quadrant  A  of  each  unit  cell  of  the 
array.  For  the  array  of  staggered  obstacles,  the  obstacles  were  placed  in  either  quadrant  A 
or  C  of  each  unit  cell,  in  accordance  to  whether  the  obstacle  row  number  (in  the  streamwise 
or  flow  direction)  was  even  or  odd,  respectively. 

The  second  obstacle  array  used  in  this  study  is  shown  in  Figure  4  and  corresponds  to  an 
aligned  array  of  31.75-mm  cubes  with  a  frontal  area  index  of  0.25  (URBOl).  The  third 
obstacle  array,  exhibited  in  Figure  5,  consists  of  31.75-mm  cubes  arranged  in  a  regular 
staggered  pattern  with  a  frontal  area  index  of  0.25  (URB02).  The  fourth  obstacle  array, 
displayed  in  Figure  6,  is  an  aligned  array  of  rectangular  blocks  with  square  cross-section  of 

velocities  in  the  streamwise  (x),  spanwise  {y)  and  vertical  (z)  directions,  respectively. 

®The  frontal  area  index  of  an  obstacle  array  is  defined  as  A/  =  Af  jA^ ,  where  Af  is  the  frontal  (windward) 
area  of  an  obstacle  and  Tl  is  the  lot  area  (surface  area  within  which  a  single  obstacle  sits  in  the  array). 
The  plan  area  index  is  defined  as  Ap  =  Ap/Al,  where  Ap  is  the  plan  (floor)  area  of  the  obstacle. 
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Figure  2:  A  schematic  diagram  of  the  geometry  of  the  Mock  Urban  Setting  Test  (MUST) 
obstacle  array. 

side  length  H  and  height  of  2H  to  give  an  array  with  a  frontal  area  index  of  0.5  (URB03). 
Figure  7  shows  the  fifth  obstacle  array  used  for  this  study,  which  corresponds  to  an  aligned 
array  of  obstacles  (cubes  or  rectangular  blocks)  with  square  cross-section  H  whose  heights 
have  been  randomly  chosen  as  H,  2H,  or  3H  (URB04).  Finally,  the  sixth  obstacle  array  is 
shown  in  Figure  8,  and  consists  of  an  aligned  array  of  rectangular  blocks  with  alternating 
rows  of  blocks  having  heights  of  2H  and  3H  resulting  in  a  frontal  area  index  of  0.625 
(URBIO). 

For  the  water  channel  simulations  involving  the  six  obstacle  arrays  summarized  above,  the 
point  source  consisted  of  a  vertical  stainless  steel  tube  (2.8  mm  LD.  and  3.1  mm  O.D.). 
The  source  emitted  a  sodium  fluorescein  dye  tracer  at  a  constant  flow  rate  of  12  x  10“^  1 
min“^  with  low  discharge  momentum  (weak  vertical  jet).  For  the  MUST  obstacle  array,  the 
point  source  was  located  at  ground  level  at  a  position  centered  in  the  street  canyon  between 
obstacles  KA  and  LA  (cf.  Figure  2).  For  arrays  URBOl  (Figure  4),  URB02  (Figure  5), 
URB03  (Figure  6)  and  URBIO  (Figure  8),  the  point  source  was  located  at  ground  level 
at  two  different  positions:  namely,  at  the  center  of  quadrant  D  behind  an  obstacle  in  a 
spanwise-oriented  street  canyon  and  at  the  center  of  quadrant  C  between  obstacles  in  a 
streamwise-oriented  street  canyon.  These  sources  were  located  in  the  unit  cell  lying  at  the 
intersection  of  the  first  row  and  eighth  column  of  obstacles  (where  the  rows  are  numbered 
in  increasing  order  in  the  streamwise  direction  from  the  leading  (windward)  edge  of  the 
array  and  the  columns  are  numbered  in  increasing  order  in  the  spanwise  direction  from  the 
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Figure  3:  Basic  unit  cell  used  to  construct  ‘urban  ’  obstacle  arrays. 


Figure  4:  Schematic  of  an  aligned  array  of  cubes  of  side  length  H  with  frontal  area  index 
Xf  =  0.25  and  plan  area  index  Xp  =  0.25  (URBOl). 

right-hand  side  of  the  array  when  looking  in  the  flow  direction  -  see  Figures  4  to  8).  For 
array  URB04  (Figure  7),  the  point  source  was  located  only  at  the  center  of  quadrant  D 
behind  an  obstacle  in  a  spanwise-oriented  street  canyon. 

The  instantaneous  concentration  field  in  the  dispersing  dye  plume  was  measured  using  the 
laser-induced  fluorescence  (LIF)  technique,  which  permitted  simultaneous  multi-point  con¬ 
centration  measurements  to  be  acquired  with  very  high  spatial  and  temporal  resolution. 
Here,  the  sodium  fluorescein  dye  was  illuminated  with  an  argon-ion  laser  beam  and  the  flu¬ 
orescent  intensities  were  measured  with  a  digital  monochrome  charge-coupled  device  (CCD) 
camera  with  a  12-bit  amplitude  resolution.  The  LIF  measurements  of  the  concentration 
were  carried  out  using  a  line-scan  configuration  with  the  line  of  laser  light  oriented  in  either 
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Figure  5:  Schematic  of  a  staggered  array  of  cubes  of  side  length  H  with  frontal  area  index 
Xf  =  0.25  and  plan  area  index  Xp  =  0.25  (URB02). 


Figure  6:  Schematic  of  an  aligned  array  of  parallelepipeds,  each  having  a  square  cross- 
section  of  side  length  H  and  a  height  of2H,  with  frontal  area  index  Xf  =  0.5  and  plan  area 
index  Xp  =  0.25  (URB03). 

a  horizontal  or  vertical  direction  through  the  dispersing  plume  to  provide  a  lateral  or  ver¬ 
tical  profile  of  the  instantaneous  plume  concentration  at  a  fixed  downstream  distance  from 
the  source  and  at  a  fixed  height  above  the  ground  surface.  Each  line  scan  consisted  of  1,016 
pixels  in  either  a  horizontal  or  vertical  line  across  the  channel,  giving  a  spatial  resolution 
between  about  0.5  mm  to  1  mm  depending  on  the  exact  positioning  of  the  camera  and  laser 
line.  The  line  scans  were  sampled  at  300  Hz.  The  sampling  times  for  each  experiment  was 
1,000  s,  yielding  N  =  300,000  individual  line  scans  per  experiment.  Each  line  scan  LIE 
measurement  provides  an  instantaneous  realization  of  the  plume  concentration  along  the 
line  of  laser  light,  and  the  light  intensity  of  each  respective  pixel  along  the  line  from  an 
entire  ensemble  of  N  line  scans  can  be  averaged  to  obtain  various  concentration  statistics 
(e.g.,  mean  concentration,  higher-order  moments  of  concentration). 

For  all  obstacle  arrays  except  URB03,  horizontal  line  scans  through  the  dispersing  array 
plume  were  made  at  a  number  of  different  downstream  distances  from  the  source,  and  for 
each  of  these  downstream  locations  a  number  of  vertical  heights  in  the  plume  were  sampled. 
For  the  MUST  array  (Figure  2),  the  horizontal  line  scans  were  made  at  hve  downstream 
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Figure  7:  Schematic  of  an  aligned  array  of  random-height  parallelepipeds,  each  having  a 
square  cross-section  of  side  length  H  and  heights  of  either  H,  2H,  or  3H  chosen  randomly, 
with  plan  area  index  Xp  =  0.25  (URB04). 


Figure  8:  Schematic  of  an  aligned  array  of  alternating-height  parallelepipeds,  each  having 
a  square  cross-section  of  side  length  H  and  heights  of  2H  or  3H  in  alternating  rows,  with 
a  frontal  area  index  Xf  =  0.625  and  plan  area  index  Xp  =  0.25  (URBIO). 

locations  (i.e.,  rows  2.5,  3.5,  4.5,  6.5,  and  9.5)  and  for  each  of  these  downstream  positions  at 
six  different  heights  (i.e.,  at  0.5,  0.75,  1.0,  1.25,  1.5,  and  2.0  times  the  obstacle  height).  Here, 
“row  2.5”  refers  to  a  streamwise  location  that  is  centered  in  the  spanwise-oriented  street 
canyon  between  rows  2  and  3  of  the  obstacle  array,  etc.  For  the  URBOl  and  URB02  arrays 
(Figures  4  and  5),  the  horizontal  line  scans  were  made  at  six  different  downstream  locations 
(rows  2.5,  3.5,  4.5,  6.5,  9.5,  and  14.5)  and  for  each  of  these  locations  at  three  different 
heights  (i.e.,  at  zjH  =  0.25,  0.5  and  1.25).  For  array  URB03,  vertical  line  scans  data 
were  acquired  for  the  dispersing  plume  along  the  mean  plume  centerline  at  six  downstream 
distances  from  the  source  (rows  2.5,  3.5,  4.5,  6.5,  9.5,  and  14.5).  Finally,  for  arrays  URB04 
and  URBIO  (Figures  7  and  8),  horizontal  line  scan  data  were  acquired  at  four  downstream 
distances  from  the  source  (rows  2.5,  3.5,  4.5,  and  9.5)  and  at  each  of  these  locations  for 
three  different  vertical  heights  (i.e.,  at  zjH  =  0.5,  1.25,  and  3.25). 
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4  Concentration  moments 


The  concentration  time  series,  extracted  from  the  measurements  of  plume  dispersion  through 
and  over  the  six  obstacle  arrays  described  in  the  previous  section,  were  processed  to  deter¬ 
mine  the  normalized  concentration  moments  up  to  order  four  (viz.,  we  calculated  {{x/C)'^) 
for  n  =  2,3,4  from  the  concentration  data),  as  well  as  the  skewness  S  and  kurtosis  K  for 
the  concentration. 

Figure  9  displays  double-logarithmic  scatterplots  (consisting  of  8,359  data  points)  of  the 
normalized  third-  and  fourth-order  concentration  moments  ((x/C)"^)  {n  =  3,4)  against  the 
normalized  second-order  concentration  moment  ((x/C)^)  and  of  kurtosis  K  against  skew¬ 
ness  S  for  the  concentration  data  measured  for  the  MUST  obstacle  array.  Figures  10  and  11 
exhibit  the  normalized  concentration  moment  diagrams  of  ((x/C)”)  (^  =  3)  4)  [consisting  of 
5,510  and  5,799  data  points,  respectively]  plotted  against  ((x/C)^)  on  a  double-logarithmic 
scale  and  of  K  plotted  against  S  for  the  aligned  array  of  cubes  (URBOl)  with  the  location  of 
the  ground-level  source  positioned  in  quadrants  C  and  D,  respectively.  Similarly,  Figures  12 
and  13  present  the  results  for  the  staggered  array  of  cubes  (URB02)  for  a  ground-level  source 
located  in  quadrants  C  and  D  consisting  of  5,538  and  5,004  data  points,  respectively;  and, 
Figures  14  and  15  show  the  normalized  concentration  moments  diagrams  for  the  aligned  ar¬ 
ray  of  parallelepipeds  of  uniform  height  2H  for  a  ground-level  source  located  in  quadrants  C 
and  D  consisting  of  938  and  1,276  data  points,  respectively.  Figure  16  presents  scatterplots 
of  the  third-  and  fourth-order  concentration  moments  versus  the  second-order  concentration 
moment  and  of  S  versus  K  for  the  aligned  array  of  random-height  parallelepipeds  (URB04) 
These  plots  consist  of  2,749  concentration  data  points.  Finally,  Figures  17  and  18  exhibit 
the  various  normalized  moment  diagrams  for  concentration  measurements  in  the  aligned  ar¬ 
ray  of  alternating-height  parallelepipeds  (URBIO)  with  the  ground-level  source  positioned 
in  quadrants  C  and  D  consisting  of  2,279  and  2,720  data  points,  respectively. 

The  points  plotted  in  Figures  9-18  cover  different  types  of  obtacles  arrays  (e.g.,  cubical 
and  non-cubical  obstacles,  different  layouts  of  obstacles  in  the  array,  arrays  with  uniform, 
random,  and  alternating  height  obstacles,  etc.)  and  a  wide  range  of  downwind,  crosswind, 
and  vertical  plume  positions  both  within  and  above  the  obstacle  array.  Even  so,  each  of 
the  scatterplots  exhibited  in  Figures  9-18  seem  to  collapse  onto  a  set  of  “universal”  curves, 
suggesting  that  this  collapse  is  a  remarkably  robust  feature  of  the  concentration  data.  This 
observed  collapse  suggests  that  there  exist  strong  correlations  between  the  various  higher- 
order  and  second-order  normalized  concentration  moments,  the  latter  of  which  also  implies 
the  existence  of  a  definitive  relationship  between  the  skewness  S  and  kurtosis  K  of  the 
concentration.  It  should  be  noted  that  there  is  more  scatter  in  the  plots  of  the  normalized 
moment  diagrams  at  the  higher  values  of  ((x/C*)^)  (corresponding  to  positions  in  the  plume 
that  are  closer  to  the  plume  edges)  and  at  higher  values  of  the  moment  order  n.  This 
increased  scatter  appears  to  be  random,  being  most  likely  attributed  to  increased  sampling 
errors  arising  from  the  measurement  of  concentration  nearer  the  plume  edges  and/or  of 
higher-order  concentration  moments. 

The  results  of  Figures  9-18  provide  compelling  evidence  that  the  concentration  PDF  for 
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(a)  {{x/Cf)  versus  {{x/Cf 


(b)  {{xlCf)  versus  ((x/C)^ 


(c)  S  versus  K 


Figure  9:  Scattergrams  of  third-  and  fourth-order  normalized  concentration  moments 
against  the  second-order  normalized  concentration  moment  and  of  skewness  S  against  kur- 
tosis  K  for  the  MUST  obstacle  array. 


a  plume  dispersing  in  a  built-up  environment  (e.g.,  arrays  of  building-like  obstacles)  can 
be  described  by  at  most  two  parameters;  namely,  a  location  and  a  scale  parameter  for 
the  PDF.  This  implies  that  one  can  use  a  simple  two-parameter  concentration  PDF  to 
predict  (approximately  or  better)  an  arbitrary  concentration  moment  given  information 
on  only  the  two  lowest  concentration  moments  [e.g.,  mean  concentration  and  mean-square 
concentration  (or,  equivalently,  concentration  variance)].  As  mentioned  earlier,  one  of  the 
goals  was  to  investigate  the  higher-order  concentration  moment  relationships  for  dispers¬ 
ing  plumes  in  order  to  suggest  an  appropriate  form  for  the  concentration  PDF  that  may 
be  used  to  formulate  a  probabilistic  model  for  the  response  of  a  chemical  agent  detector. 
Towards  this  objective,  we  consider  two  simple  concentration  PDF  models  (whose  form  is 
completely  specified  by  two  parameters)  and  examine  their  flexibility  in  representing  the 
observed  higher-order  concentration  moment  relationships  measured  in  a  built-up  environ¬ 
ment;  namely,  the  clipped-gamma  and  intermittent  exponential  PDFs  [cf.  Eqs.  (2)  and  (15), 
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s 

(c)  S  versus  K 

Figure  10:  Scatter  grams  of  third-  and  fourth-order  normalized  concentration  moments 
against  the  second-order  normalized  concentration  moment  and  of  skewness  S  against 
kurtosis  K  for  the  URBOl  obstacle  array  with  the  source  located  in  quadrant  C  (in  a 
streamwise-oriented  street  canyon  between  two  columns  of  obstacles). 


respectively] . 

For  comparative  purposes,  the  theoretical  relationships  between  ((x/C)'^)  (n  =  3,4)  and 
((x/C*)^)  and  between  S  and  K,  generated  by  the  clipped-gamma  PDF  and  intermittent 
exponential  PDF,  have  been  superimposed  on  the  scatterplots  in  Figures  9-18.  In  all 
these  figures,  the  solid  and  dashed  curves  represent  the  theoretical  curves  predicted  using 
the  clipped-gamma  and  intermittent  exponential  distributions,  respectively.  It  is  seen  that 
both  the  intermittent  exponential  and  clipped-gamma  distributions  appear  to  provide  rather 
good  fits  to  the  measurements  over  the  full  range  of  conditions  covered  by  the  concentration 
data  sets.  Even  so,  it  appears  that  generally  the  concentration  moments  are  slightly  under¬ 
predicted  by  the  intermittent  exponential  distribution  (viz.,  more  of  the  data  points  in  the 
concentration  moment  diagrams,  and  especially  in  the  plots  of  kurtosis  against  skewness. 
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(c)  S  versus  K 

Figure  11:  Scatter  grams  of  third-  and  fourth-order  normalized  concentration  moments 
against  the  second-order  normalized  concentration  moment  and  of  skewness  S  against  kur- 
tosis  K  for  the  URBOl  obstacle  array  with  the  source  located  in  quadrant  D  (in  a  spanwise- 
oriented  street  canyon  between  two  rows  of  obstacles). 


lie  slightly  above  the  curve  generated  by  the  exponential  PDF). 

The  observed  concentration  moment  relationships  for  plumes  dispersing  in  various  obstacle 
arrays  (summarised  in  Figures  9-18)  are  in  good  agreement  with  predictions  provided  by 
the  clipped-gamma  PDF.  Interestingly,  Yee  and  Chan  [14]  demonstrated  that  the  clipped- 
gamma  PDF  provided  predictions  of  concentration  moment  relationships  that  were  in  very 
good  conformance  with  the  measured  concentration  moment  relationships  obtained  from 
the  extensive  CONFLUX  experiments.  These  experiments  provided  concentration  data  for 
plumes  dispersing  over  level,  unobstructed  terrain  that  covered  a  very  wide  range  in  plume 
positions  and  atmospheric  conditions.  Taken  together,  this  implies  that  the  concentration 
moment  relationships  for  plumes  dispersing  in  the  complex  urban  terrain  are  identical  to 
those  for  plumes  dispersing  in  simple  (rural)  terrain. 
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(c)  S  versus  K 

Figure  12:  Scattergrams  of  third-  and  fourth-order  normalized  concentration  moments 
against  the  second-order  normalized  concentration  moment  and  of  skewness  S  against 
kurtosis  K  for  the  URB02  obstacle  array  with  the  source  located  in  quadrant  C  (in  a 
streamwise-oriented  street  canyon  between  two  columns  of  obstacles). 


It  is  known  that  large  groups  of  obstacles  (characteristic  of  a  built-up  environment)  can 
have  a  profound  effect  on  the  mean  concentration  and  concentration  variance  in  a  dispersing 
plume  in  comparison  to  plume  dispersion  over  open  terrain.  In  an  obstacle  array,  non- 
Gaussian  distributions  of  mean  concentration  are  possible  due  to  the  splitting  of  plumes  into 
two  by  the  horseshoe  vortex  that  may  form  around  the  base  of  an  obstacle.  There  may  be  an 
effective  lifting  of  the  mean  plume  centroid  caused  by  deflection  of  mean  flow  streamlines 
over  an  array  of  obstacles.  The  trajectory  of  the  mean  plume  centerline  in  the  lateral 
direction  may  be  offset  from  the  direction  of  the  mean  wind  aloft  due  to  crosswind  channeling 
of  the  flow  in  the  array  of  obstacles.  In  certain  circumstances,  channeling  of  the  flow  along 
array  passages  may  inhibit  the  lateral  spread  of  the  mean  plume.  Alternatively,  topological 
diffusion  of  the  plume  due  to  streamline  divergence  can  lead  to  an  enhanced  lateral  mean 
plume  spread,  whereas  rapid  vertical  mixing  in  the  recirculating  wakes  behind  obstacles 
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(c)  S  versus  K 

Figure  13:  Scattergrams  of  third-  and  fourth-order  normalized  concentration  moments 
against  the  second-order  normalized  concentration  moment  and  of  skewness  S  against  kur- 
tosis  K  for  the  URB02  obstacle  array  with  the  source  located  in  quadrant  D  (in  a  spanwise- 
oriented  street  canyon  between  two  rows  of  obstacles). 

may  lead  to  an  enhanced  vertical  mean  plume  spread.  These  important  modifications  of 
the  structure  of  the  mean  concentration  field,  in  a  plume  dispersing  through  an  array  of 
obstacles,  have  been  discussed  by  Davidson  et  al  [23],  Macdonald  et  al  [24],  and  Yee  and 
Biltoft  [17]. 

Yee  and  Biltoft  [17]  demonstrated  that  the  second  moment  of  concentration  (or,  concen¬ 
tration  variance)  in  a  plume  dispersing  through  an  array  of  obstacles  is  reduced,  relative 
to  that  observed  in  a  comparable  plume  dispersing  over  open  terrain.  This  reduction  in 
the  concentration  variance  is  the  consequence  of  two  physical  mechanisms:  (1)  a  significant 
reduction  in  the  meandering  of  the  instantaneous  plume  in  an  obstacle  array  (relative  to 
an  open-terrain  plume)  due  to  an  observed  increase  in  the  lateral  spreads  of  the  plume  and 
a  reduction  in  the  scale  of  turbulence  between  the  obstacles  in  the  array,  and  (2)  a  signif¬ 
icant  reduction  in  concentration  variance  of  the  in-plume  fluctuations  in  the  array  plume 
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s 

(c)  S  versus  K 

Figure  14:  Scattergrams  of  third-  and  fourth-order  normalized  concentration  moments 
against  the  second-order  normalized  concentration  moment  and  of  skewness  S  against 
kurtosis  K  for  the  URB03  obstacle  array  with  the  source  located  in  quadrant  C  (in  a 
streamwise-oriented  street  canyon  between  two  columns  of  obstacles). 


(in  comparison  to  the  open-terrain  plume)  owing  to  the  enhanced  small-scale  mixing  of  the 
plume  material  in  the  high  intensity  turbulence  that  is  characteristic  of  the  flow  within  an 
obstacle  array. 

In  spite  of  the  significant  modifications  in  the  structure  of  a  plume  dispersing  in  an  obstacle 
array  when  compared  to  an  open-terrain  plume,  the  modifications  manifest  themselves  in 
their  effects  on  the  first-  and  second-order  concentration  moments  (or,  equivalently,  on 
the  mean  concentration  and  concentration  variance).  Interestingly,  Figures  9-18  seem  to 
suggest  that  despite  the  fact  that  an  obstacle  array  has  a  significant  effect  on  the  first  two 
moments  of  the  concentration,  the  relationships  of  the  third-  and  higher-order  normalized 
concentration  moments  to  the  second-order  normalized  concentration  moment  in  an  array 
plume  are  exactly  the  same  as  those  observed  in  an  open-terrain  plume.  This  suggests  that 
the  plume  concentration  fluctuation  structure  in  an  obstacle  array  is  the  same  as  that  which 
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s 

(c)  S  versus  K 

Figure  15:  Scattergrams  of  third-  and  fourth-order  normalized  concentration  moments 
against  the  second-order  normalized  concentration  moment  and  of  skewness  S  against  kur- 
tosis  K  for  the  URB03  obstacle  array  with  the  source  located  in  quadrant  D  (in  a  spanwise- 
oriented  street  canyon  between  two  rows  of  obstacles). 


would  be  observed  in  an  open-terrain  plume  at  a  much  greater  distance  downwind  when  the 
instantaneous  plume  has  grown  to  fill  the  mean-plume  width  (hence,  dramatically  reducing 
the  contribution  of  large-scale  meander  to  the  concentration  fluctuations)  and  when  the 
continuous  molecular  mixing  has  had  a  chance  to  smooth  out  the  in-plume  fluctuations. 
The  latter  effect  on  the  structure  of  an  open-terrain  plume  is  identical  to  the  effect  of  the 
more  rapid,  homogeneous  mixing  in  an  array  plume  due  to  the  small-scale,  high-intensity 
turbulence  that  is  characteristic  of  the  flow  within  an  obstacle  array.  This  explains  why  the 
concentration  moment  relationships  observed  in  an  array  plume  are  identical  (approximately 
or  better)  to  those  observed  in  an  open-terrain  plume. 
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(b)  {{x/Cf)  versus  {{x/CY) 


s 

(c)  S  versus  K 


Figure  16:  Scatter  grams  of  third-  and  fourth-order  normalized  concentration  moments 
against  the  second-order  normalized  concentration  moment  and  of  skewness  S  against  kur- 
tosis  K  for  the  URB04  obstacle  array  with  the  source  located  in  quadrant  D  (in  a  spanwise- 
oriented  street  canyon  between  two  rows  of  obstacles). 


5  Concentration  PDF 


In  this  section,  we  compare  the  shape  of  the  model  probability  distributions  to  the  observed 
concentration  data  at  a  number  of  plume  locations  for  two  of  the  obstacle  arrays.  The 
concentration  probability  distribution  embodies  the  information  required  to  predict  the 
probability  of  exposure  to  high  concentrations  and,  as  such,  contains  the  information  on 
all  the  moments  of  the  fluctuating  plume  concentration.  This  information  can  be  used  to 
derive  all  one-point  statistics  of  the  concentration  field. 

To  this  end,  we  extract  the  cumulative  and  exceedance  probability  distributions  from  the 
array  plume  concentration  data  and  compare  these  measured  distributions  to  two  model 
cumulative  and  exceedance  probability  distributions;  namely,  to  the  clipped-gamma  and 
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intermittent  exponential  distributions.  In  view  of  Eqs.  (2)  and  (15),  the  cumulative  distri¬ 
bution  function  for  the  concentration  at  a  fixed  receptor  point  can  be  written  as 

F{c;x,t)  =  7F+(c;x,t)  +  {1  -  7)F(c;x,t),  (22) 

where  F~^(c]x,t)  is  the  probability  distribution  for  the  in-plume  concentration  fluctuations 
and  H{c]x,t)  is  a  unit  step  function  with  a  step  at  c  =  0  of  magnitude  (1  —  7),  which 
expresses  the  probability  of  observing  zero  concentration  at  receptor  location  x  at  time  t. 
For  the  clipped-gamma  distribution  [cf.  Eq.  (2)], 

F(c;x,t)  = 


Similarly,  for  the  intermittent  exponential 

F{c;:s.,t)  = 

where  7  is  determined  by  Eq.  (19).  Observe  from  Eqs.  (23)  and  (24)  that 

F(0";x,t)  =  1  -  7,  (25) 

is  the  probability  of  observing  zero  concentration  at  receptor  location  x  and  time  t. 

Figure  19  presents  measured  cumulative  distribution  functions  (CDFs)  of  the  normalized 
concentration  c/C  obtained  at  five  different  downstream  locations  along  the  mean  plume 
centerline  at  half  canopy  height  for  dispersion  in  the  MUST  obstacle  array.  The  clipped- 
gamma  (solid  line)  and  intermittent  exponential  (dashed  line)  distributions,  with  the  same 
normalized  mean-squared  concentration  ((x/C*)^)  as  the  measured  result,  are  also  exhibited 
in  Figure  19  for  comparison.  The  clipped-gamma  distribution  is  seen  to  be  to  in  excellent 
conformance  with  the  measured  concentration  CDFs  at  the  various  downstream  locations. 
However,  the  intermittent  exponential  distribution  provides  a  very  poor  fit  to  the  mea¬ 
sured  concentration  CDFs.  More  specifically,  the  lower  tail  of  the  intermittent  exponential 
distribution  is  much  heavier  than  that  for  the  observed  concentration  CDF  (viz.,  the  inter¬ 
mittent  exponential  distribution  predicts  a  larger  probability  for  non-exceedances  of  small 
concentration  levels  with  c/C  <  1  than  that  given  by  the  observed  values). 

Generally,  for  a  number  of  practical  applications,  it  is  the  largest  peak  concentrations 
that  are  of  greatest  interest.  Because  it  is  the  prediction  of  the  likelihood  of  extreme 
events  that  is  important  in  the  hazard  assessment  of  toxic  gas  releases,  it  is  important  to 
examine  the  upper  tail  of  the  concentration  distributions.  To  that  purpose,  we  display  in 
Figure  20  the  exceedance  probability  distributions,  1  —  F{c/C),  for  the  same  plume  locations 
shown  in  Figure  19  for  dispersion  in  the  MUST  obstacle  array.  The  exceedance  probability 
distributions  have  been  plotted  on  a  logarithmic  scale  in  order  to  emphasize  the  upper  tails. 
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(a)  {ix/Cf)  versus  {{x/Cf)  (b)  {{x/Cf)  versus  {{x/CY) 


s 

(c)  S  versus  K 


Figure  17:  Scattergrams  of  third-  and  fourth-order  normalized  concentration  moments 
against  the  second-order  normalized  concentration  moment  and  of  skewness  S  against 
kurtosis  K  for  the  URBIO  obstacle  array  with  the  source  located  in  quadrant  C  (in  a 
streamwise-oriented  street  canyon  between  two  columns  of  obstacles). 


Figure  20  indicates  that  the  clipped-gamma  distribution  predicts  the  upper  tail  extremely 
well.  In  contrast,  the  intermittent  exponential  distribution  exhibits  considerably  longer 
upper  tails  than  the  observed  concentration  probability  distribution.  This  implies  that  the 
intermittent  exponential  distribution  tends  to  overpredict  (significantly)  the  probability  of 
occurrence  of  large  concentrations  (c/C  >  1). 

As  another  example,  Figures  21  and  22  exhibit  the  measured  CDFs  and  EDFs,  respec¬ 
tively,  of  the  normalized  concentration  c/C  obtained  at  six  different  downstream  locations 
at  half  canopy  height  (z/H  =  0.5)  along  the  mean  centerline  for  the  plume  dispersing  in  the 
URBOl  obstacle  array  with  the  source  located  in  quadrant  D.  Additionally,  these  figures 
show  the  curves  provided  by  the  clipped-gamma  (solid  line)  and  intermittent  exponential 
(dashed  lines)  distributions,  which  were  generated  using  the  measured  mean-square  concen¬ 
tration  {(x/CY)  obtained  at  each  of  these  plume  locations.  Again,  it  can  be  seen  that  the 
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(a)  {{x/Cf)  versus  {{x/Cf 


(b)  {{xlCf)  versus  ((x/C)^ 


s 


(c)  S  versus  K 


Figure  18:  Scattergrams  of  third-  and  fourth-order  normalized  concentration  moments 
against  the  second-order  normalized  concentration  moment  and  of  skewness  S  against  kur- 
tosis  K  for  the  URBIO  obstacle  array  with  the  source  located  in  quadrant  D  (in  a  spanwise- 
oriented  street  canyon  between  two  rows  of  obstacles). 


clipped-gamma  distribution  provides  very  good  predictions  for  the  shape  of  the  measured 
concentration  probability  distributions.  More  importantly,  the  extent  of  the  upper  tails  of 
the  concentration  probability  distributions  are  predicted  very  well  using  the  clipped-gamma 
distribution.  In  contrast,  the  intermittent  exponential  distribution  is  seen  to  provide  rather 
poor  predictions  for  the  measured  concentration  probability  distribution.  This  model  dis¬ 
tribution  tends  to  overpredict  the  probability  of  non-exceedances  of  small  concentration 
levels  (c/C  <  1)  and  exceedances  of  large  concentration  levels  (c/C  >  1). 
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F{c/C)  F(c/C) 


c/C 


(a)  Plume  centerline  at  row  2.5 


c/C 


(c)  Plume  centerline  at  row  4.5 


c/C 


(b)  Plume  centerline  at  row  3.5 


c/C 


(d)  Plume  centerline  at  row  6.5 


c/C 


(e)  Plume  centerline  at  row  9.5 

Figure  19:  Cumulative  distribution  function  (CDF),  F{c/C),  of  the  normalized  concentra¬ 
tion  measured  at  various  receptor  positions  at  half  canopy  height  along  the  mean  centerline 
of  a  plume  dispersing  in  the  MUST  obstacle  array. 


26 


DRDC  Suffield  TR  2008-077 


Figure  20:  Exceedance  distribution  function  (EDF),  1  —  F{c/C),  of  the  normalized  con¬ 
centration  measured  at  various  receptor  positions  at  half  canopy  height  along  the  mean 
centerline  of  a  plume  dispersing  in  the  MUST  obstacle  array. 
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F(c/0  F(c/C)  F{c/C) 


c/C 

(b)  Plume  centerline  at  row  3.5 


c/C 


c/C 


(c)  Plume  centerline  at  row  4.5 


(d)  Plume  centerline  at  row  6.5 


c/C 

(e)  Plume  centerline  at  row  9.5 


Figure  21 :  Cumulative  distribution  function  (CDF),  F{c/C),  of  the  normalized  concentra¬ 
tion  measured  at  various  receptor  positions  at  half  canopy  height  along  the  mean  centerline 
of  a  plume  dispersing  in  the  URBOl  obstacle  array  with  the  source  in  quadrant  D. 
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Figure  22:  Exceedance  distribution  function  (EDF),  1  —  F{c/C),  of  the  normalized  con¬ 
centration  measured  at  various  receptor  positions  at  half  canopy  height  along  the  mean 
centerline  of  a  plume  dispersing  in  the  URBOl  obstacle  array  with  the  source  in  quadrant 
D. 


DRDC  Suffield  TR  2008-077 


29 


6  Implications  for  sensor  modeiing 


A  critical  capability  gap  in  current  emergency  and  retrospective  management  efforts,  di¬ 
rected  at  terrorist  incidents  involving  the  covert  release  of  a  chemical,  biological  and  radio¬ 
logical  (CBR)  agent  in  a  densely  populated  urban  center,  is  the  localization  of  the  unknown 
source(s)  following  event  detection  by  an  array  of  independent  CBR  sensors.  A  consider¬ 
able  research  effort  has  been  undertaken  in  the  past  two  years  to  develop  a  methodology 
to  address  the  inverse  prediction  of  agent  source  location  and  strength  using  a  finite  num¬ 
ber  of  noisy  concentration  measurements  obtained  from  an  array  of  CBR  detectors.  To 
this  purpose,  Bayesian  probabilistic  inference  has  been  applied  successfully  to  the  problem 
of  source  resconstruction  in  complex  environmental  conditions  (e.g.,  complex  and/or  ur¬ 
ban  terrain  involving  highly  disturbed  wind  fields  exhibiting  a  significant  degree  of  spatial 
inhomogeneity  and/or  temporal  non-stationarity)  [25],  [26],  [27]. 

In  the  studies  of  source  reconstruction  cited  here,  the  time-averaged  (or,  ensemble)  mean 
concentration  measured  at  the  sensor  location  was  used  for  the  source  reconstruction,  rep¬ 
resenting  the  mean  of  either  a  number  of  realizations  of  concentration  or  concentrations 
averaged  over  periods  of  time  of  several  minutes  or  hours.  However,  many  detectors  pro¬ 
vide  measurements  of  concentration  occurring  on  time  scales  of  a  few  seconds  or  less,  not 
the  average  over  minutes  or  hours.  Moreover,  many  of  these  detectors  only  provide  short 
time-averaged  concentrations  that  have  been  discretized  into  a  small  number  of  amplitude 
levels  (or,  bars).  For  example,  the  Chemical  Agent  Monitor  (CAM)  which  uses  ion  mobility 
spectrometry  (Smiths  Detection)  and  the  AP4C  chemical  detector  which  uses  flame  spec¬ 
trophotometry  (Proengin  Inc.)  are  hand-held  chemical  warfare  agent  field  detectors  that 
report  short  time-averaged  concentration  levels  discretized  on  an  eight-  and  five-bar  scale, 
respectively. 

In  order  to  use  this  type  of  concentration  information  for  source  reconstruction,  we  need  a 
description  of  the  concentration  fluctuations,  as  well  as  the  ensemble  mean  concentration 
in  the  formulation  of  the  problem.  Indeed,  a  chemical  agent  detector  with  a  short  aver¬ 
aging  time  (less  than  tens  of  seconds)  placed  at  a  fixed  position  downwind  of  a  release  of 
a  contaminant  will  experience  concentration  fluctuations  because  of  plume  sweeping  and 
small-scale  turbulent  mixing.  These  fluctuations  can  deviate  significantly  from  the  mean 
concentration  value,  and  models  for  the  response  of  these  sensors  must  account  properly 
for  the  effects  of  random  concentration  fluctuations  in  hazardous  gas  releases.  The  purpose 
of  this  section  is  to  show  how  this  can  be  accomplished. 

A  key  quantity  required  for  Bayesian  source  inversion  is  the  likelihood  function  which  en¬ 
codes  all  the  information  provided  by  the  concentration  data  about  the  unknown  contam¬ 
inant  source.  The  likelihood  function  provides  a  measure  of  how  likely  the  concentration 
measurements  made  by  a  detector  (e.g.,  CAM,  AP4C)  would  be  for  a  given  source  dis¬ 
tribution.  We  will  now  show  how  to  formulate  a  likelihood  function  for  a  chemical  agent 
detector,  such  as  the  CAM  or  AP4C,  in  terms  of  the  concentration  PDF  which  provides  the 
one-point  statistical  information  about  the  instantaneous  plume  or  cloud  concentration. 

A  chemical  agent  bar  sensor  will  exhibit  a  bar  level  bi  (i  =  1,2, . . .  ,  Nh  where  Nf,  is  the 
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maximum  number  of  bar  levels),  if  the  measured  concentration  at  the  receptor  location  x 
averaged  over  a  short  time  window  (typically  between  about  5-10  s)  centered  at  time  t,  lies 
in  the  concentration  semi-open  interval  Here,  cf  and  cf  are  the  lower  and  upper 

bounds,  respectively,  of  the  concentration  that  defines  the  i-th  bar  level  for  a  chemical  agent 
detector.  In  other  words,  the  chemical  agent  detector  will  issue  a  “reading”  for  the  i-th  bar 
level  bi  if  the  measured  concentration  Cm  G  (cf ,  cf].  Now,  we  can  apply  the  product  rule  of 
probability  calculus  to  obtain  the  following  expression  for  the  probability  that  a  bar  level 
bi  is  observed  by  a  chemical  agent  detector,  conditioned  on  a  release  from  a  given  source 
distribution  5:® 

roo  poo  poo  _  _  _ 

V{bi)=  /  /  /  /  f{cm\c)f{c\C,c^)f{C\S)f{c^\S)dc^dCdcdcm,  (26) 

Jcf  Jo  Jo  Jo 

for  i  =  1, 2, . . . ,  Nh.  In  Eq.  (26),  f{cm\c)  is  the  probability  density  function  of  the  measured 
concentration  Cm  conditioned  on  the  true  (unknown)  concentration  c;  f{c\C,c^)  is  the  con¬ 
centration  PDF,  which  is  implicitly  assumed  to  depend  on  two  parameters  (namely,  the 
mean  concentration  C  and  the  mean-square  concentration  f{C\S)  is  the  PDF  for  C 
given  the  source  distribution  S;  and,  f{c^\S)  is  the  PDF  for  given  the  source  distribution 
S.  Next,  we  consider  the  specification  or  assignment  of  each  of  these  component  PDFs. 

First,  let  us  consider  the  assignment  of  f{cm\c).  If  the  true  (unknown)  concentration  is  c, 
then  the  measured  concentration  datum  Cm  is  given  by 

Cjn  =  c  +  ei,  (27) 

where  ei  represents  the  measurement  (detector)  noise,  which  includes  the  effects  of  digitiza¬ 
tion  errors,  calibration  errors,  and  shot  noise.  It  is  assumed  that  the  measurement  error  ei 
is  known  to  have  zero  mean  and  variance  af .  The  principle  of  maximum  entropy  provides  a 
powerful  criterion  for  formulation  of  a  prior  noise  probability  density  function  that  encodes 
our  state  of  knowledge  about  the  measurement  noise  ei  (see  Jaynes  [28]  for  an  insightful 
discussion  of  the  maximum  entropy  principle).  Application  of  the  principle  of  maximum 
entropy  informs  us  that  a  Gaussian  distribution  is  the  most  conservative  choice  for  the  PDF 
of  ei,  or  equivalently,  for  /(cm|c); 

/(Cm|c)  =  f{cm\c,  cjf)  =  exp  ^  ,  (28) 

where  the  explicit  dependence  of  f{cm\c)  on  af  has  been  indicated.  It  should  be  noted  that 
any  other  noise  PDF  will  have  more  compact  support  (lower  entropy)  than  the  Gaussian 
given  the  finite  noise  power  af;  thus,  the  assignment  of  the  Gaussian  form  for  the  noise 
PDF  will  give  the  most  conservative  representation  for  the  structure  of  the  measurement 
noise. 

®The  dependence  of  various  quantities  on  x  and  t  has  been  suppressed  to  simplify  the  notation.  Fur¬ 
thermore,  the  conditioning  of  the  probability  p{bi)  on  the  source  distribution  S  has  been  suppressed  in  the 
notation. 

^This  corresponds  to  a  slight  abuse  of  the  notation  used  earlier,  in  which  the  mean-square  concentration 
was  denoted  as  [cf.  Eq.  (7)].  Hence,  in  this  section,  the  identification  =  (y^)  should  be  made. 
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Next,  let  us  focus  on  specification  of  the  concentration  PDF  f{c\C,c^).  In  a  previous 
section,  it  was  shown  that  the  form  of  the  concentration  PDF  can  be  specified  using  only 
two  parameters,  and  that  the  clipped-gamma  distribution  fits  the  concentration  observations 
very  well.  In  consequence,  we  can  assign  a  clipped-gamma  distribution  for  f{c\C,c^).  By 
virtue  of  Eqs.  (2),  (3)  and  (10),  the  concentration  PDF  f{c\C,  c^)  can  be  explicitly  expressed 
as  follows:® 


/(c|C,c2) 


c-FA\^  ^  exp(-(c A)/s) 
s  )  sT{k) 


+  (1  -7)'5(c), 


(29) 


where 


7  =  mm  (1,  -=^ 


r(fc;A/s) 

m 


(30) 


It  should  be  stressed  that  k,  s,  A  is  uniquely  determined  from  a  knowledge  of  C  and  c^. 
More  specihcally,  the  latter  quantities  determine  the  normalized  mean-square  concentration 
(?  =  ((x/C)^))  and  as  exhibited  in  Figure  1,  the  parameters  k,  s,  and  A  of  the  clipped- 

gamma  distribution  are  uniquely  specified  given  the  value  of  (?  . 


Next,  let  us  consider  f{C\S).  The  true  mean  concentration  C,  at  the  location  of  the  chemical 
agent  detector,  resulting  from  the  release  of  a  contaminant  from  a  given  source  distribution 
S  is  not  known.  This  quantity  needs  to  be  estimated  using  an  atmospheric  dispersion  model 
(e.g.,  C  can  be  estimated  using  either  a  Gaussian  puff  or  a  one-particle  Lagrangian  stochastic 
trajectory  simulation  model).  Hence,  f{C\S)  specifies  the  uncertainty  in  the  estimation  of 
C  using  a  dispersion  model  for  a  hazardous  gas  release  from  a  given  source  distribution 
S.  To  this  purpose,  we  first  assume  that  the  predicted  (or,  model)  mean  concentration  Cp, 
obtained  from  an  atmospheric  dispersion  model,  provides  an  unbiased  estimate  for  the  true 
(unknown)  mean  concentration  C.  The  one  standard  deviation  uncertainty  in  the  estimate 
Cp  for  C  is  iT2,  the  latter  of  which  is  assumed  to  be  known.®  This  uncertainty  arises  from 
model  errors  in  the  representation  of  various  physical  processes  in  the  dispersion  model  used 
to  predict  the  mean  concentration  and  input  errors  in  the  values  of  empirical  parameters 
and/or  specification  of  the  input  meteorology  (initial  and  boundary  conditions)  used  by  the 
dispersion  model.  Next,  we  assume  that  f{C\S)  can  be  specified  by  a  gamma  distribution 
with  mean  Cp  and  standard  deviation  <72,  so 


f{C\S)  =  f{C\S,Cp,a2) 


exp(-C'C'p/cj|) 

{allCpf'l^hiClIal)  ’ 


(31) 


where  the  explicit  dependence  of  f{C\S)  on  Cp  and  cj2  has  been  indicated. 


Finally,  we  consider  the  assignment  of  f{c^\S)  which  encodes  the  uncertainty  in  our  predic¬ 
tion  of  the  mean-square  concentration  c^p  for  a  hazardous  gas  release  from  a  given  source 
distribution  S,  the  true  value  of  the  mean-square  concentration  being  unknown.  The 
mean-square  concentration  can  be  modeled  using  a  fluctuating  plume  model  [29],  [30]  or 

^Recall  that  for  the  clipped-gamma  distribution,  when  c^/D^  <  3  the  intermittency  factor  is  set  equal  to 
unity  identically  and  this  distribution  reduces  to  the  usual  gamma  distribution  under  these  circumstances. 

^Normally,  the  expected  uncertainty  in  the  prediction  of  the  mean  concentration  using  a  dispersion  model 
is  available  from  model  validation  (evaluation)  exercises  conducted  with  benchmark  dispersion  data  sets. 
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a  second-order  closure  model  [31],  [21],  The  predicted  (model)  mean-square  concentration 
c^p  is  assumed  to  be  an  unbiased  estimate  for  the  true  (unknown)  mean-square  concen¬ 
tration  c^.  The  one  standard  deviation  uncertainty  in  the  estimate  c^p  for  c?  is  CJ3,  the 
latter  providing  the  scale  of  the  uncertainty  in  the  determination  of  the  true  mean-square 
concentration  arising  from  model  and  input  errors.  Furthermore,  it  is  assumed  that  /(c^jS) 
can  be  specified  by  gamma  distribution  with  mean  c^p  and  standard  deviation  1T3,  so 


f{c^\S)  =  /(c2|5,c2p,cr3) 


— _ t  2,2 


exp(— 


/^i) 


(32) 


where  the  explicit  dependence  of  f{c^\S)  on  c^p  and  CJ3  has  been  indicated. 


Now,  we  can  insert  Eq.  (28)  into  Eq.  (26)  and  integrate  with  respect  to  Cm  to  simplify  the 
form  for  p{bi).  This  gives 


/(c|C,c2)/(C|5)/(c2|5)  dc2  dCdc, 


(33) 


where  erf(x)  denotes  the  error  function.  Usually,  the  model  and  input  errors  in  the  pre¬ 
diction  of  C  and  are  much  larger  than  the  errors  in  the  measurement  of  c.  In  this  case, 
a  I  <C  CJ2,  CJ3,  so  that  f{cm\c)  in  Eq.  (28)  can  be  considered  to  be  so  sharply  peaked  about  the 
true  concentration  c  that  the  function  can  be  reasonably  approximated  as  a  delta  function: 

/(cmjc)  =  lim  exp  =d{cm-c),  (34) 

cri^0+  y'2TTai  V  / 

so  that  as  fii  ^  O'*'  in  Eq.  (28),  the  approximation  of  /(cmjc)  using  a  delta  function  becomes 
exact.  With  the  approximation  for  f{cm\c)  given  by  Eq.  (34),  the  form  for  p{bi)  given  by 
Eq.  (26)  simplifies  to  (on  noting  that  Cm  ^  c  for  this  case) 

poo  poo  _  _  _ 

P{bi)=  I  /  /(c|C,?)/(C|5)/(?|5)d?dCdc.  (35) 

Jc<  Jo  Jo 


Eor  the  case  where  the  measurement  error  is  assumed  to  be  neglible  with  respect  to  the 
model  and  input  errors,  insertion  of  the  functional  form  for  the  concentration  PDE  given 
by  Eq.  (29)  into  Eq.  (35)  and  performing  the  integration  with  respect  to  c  yields 


P{bi) 


pOiO  poo 

'T{k-ic<+X)/s)-T{k-,ic>  +  X)/sy 

Jo  Jo 

m 

fiC\S)fic^\S)  dd^  dC, 


(36) 


for  i  =  1,  2, . . . ,  Nf,.  The  remaining  integrations  with  respect  to  C  and  c2  in  Eq.  (36)  cannot 
be  evaluated  analytically  and  will  need  to  be  evaluated  numerically.  Note  in  Eq.  (36)  that 
k,  s,  and  A  are  implicit  functions  of  C  and  c2  (more  specifically,  of  ((x/C*)^)  =  (?  jC’^')  as 
shown  in  Eigure  1. 
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The  bar  levels  6*  (i  =  1,2,  ...jA/fe)  correspond  to  non-zero  bar  “readings”.  A  zero  bar 
reading,  which  we  designate  using  bo,  corresponds  to  the  case  when  the  concentration 
level  is  less  than  cf  (lower  bound  defining  the  first  bar  level).  For  consistency  with  our 
previous  notation  for  non-zero  bar  levels,  the  zero  bar  “reading”  would  be  associated  with 
the  probability  that  the  concentration  c  lies  in  the  interval  Cq  ^  0“  and  Cq  cf .  For  this 
special  case,  the  Dirac  delta  function  term  in  Eq.  (29)  would  contribute  to  the  integral  over 
c  in  Eq.  (35).  Consequently,  for  this  case,  the  probability  that  the  detector  would  report  a 
zero  bar  reading  (corresponding  to  a  “nominal”  zero  concentration)  is  given  by 

pM  =  -r(Mc<  +  A)/.)j 

roo  roo  _  _ 

+  /  /  {l-^)f{C\S)f(^\S)(0dC,  (37) 

Jo  Jo 

where  7  is  defined  by  Eq.  (30).  The  double  integrations  in  Eq.  (37)  will  need  to  be  evaluated 
numerically. 

It  should  be  noted  that  p{bi)  {i  =  0,1,2,...  ,N(i)  given  by  Eqs.  (36)  and  (37)  should  be 
identified  with  the  likelihood  C{bi\S)  that  a  detector  reports  the  bar  level  6*  for  a  hazardous 
gas  release  associated  with  a  given  source  distribution  S.  This  likelihood  function  accounts 
fully  for  uncertainties  in  model  predictions  for  C  and  c^,  required  in  the  specification  of  the 
concentration  PDE.  More  specifically  and  succinctly,  C{bi\S)  =  p{bi)  {i  =  0, 1,2, . . .  ,Nb). 
This  quantity  can  be  used  for  the  likelihood  function  in  a  Bayesian  source  inversion  scheme 
involving  concentration  measurements  made  by  a  chemical  agent  detector  that  outputs  only 
a  discrete  level  bar  reading  (rather  than  a  continuous  level  concentration  reading). 

It  was  not  mentioned  previously,  but  the  determination  of  p{bi)  given  by  Eq.  (26)  requires 
information  on  the  concentration  moments  C  and  as  perceived  by  the  detector  (viz., 
the  concentration  moments  as  measured  by  a  detector  with  the  averaging  time  Ta)-  Time¬ 
averaging,  such  as  that  imposed  by  a  slower  response  detector,  does  not  affect  the  mean 
concentration  C  (viz.,  C  is  invariant  to  a  time-averaging  operation).  However,  the  time¬ 
averaging  inherent  in  the  detector  response  will  affect  the  mean-square  concentration  c^, 
owing  to  the  fact  that  time-averaging  will  smooth  out  the  fine  structure  of  rapid  fluctuations 
in  the  concentration  time  series  measured  by  the  detector. 

The  effects  of  averaging  time  on  the  expected  reduction  in  the  concentration  variance 
fj2  =  q2  _  (j2  been  discussed  in  Hanna  [32]  —  namely,  the  variance  of  concentra¬ 
tion  fluctuations  at  a  receptor  point,  averaged  over  time  Tq,  is  related  to  the  concentration 
variance  with  no  time  averaging  (viz.,  zero  averaging  time)  as  follows: 

=  2—  1  -  —(1  -  exp(-Ta/rc))  ,  (38) 

where  is  the  concentration  variance  perceived  by  a  detector  with  an  intrinsic  averaging 
time  Ta-  In  Eq.  (38),  Tc  is  the  integral  time  scale  of  the  turbulent  concentration  fluctuations 
defined  as 

poo 

Tc=  R{t)  (It,  (39) 

Jo 
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where  R{t)  is  the  autocorrelation  function  of  the  concentration  fluctuations  (r  is  the  cor¬ 
relation  lag  time). 

A  parameterization  of  the  integral  scale  for  concentration  fluctuations  has  been  suggested 
by  Sykes  [33].  The  effects  of  time-averaging  imposed  by  a  detector  can  be  accounted  for 
by  using  a  model  to  predict  the  unaveraged  concentration  variance  (or,  equivalently, 
the  unaveraged  mean-square  concentration  c^p),  and  then  using  Eq.  (38)  to  compute  the 
concentration  variance,  c'^plra)  appropriate  for  the  averaging  time  Ta  imposed  by  the  detec¬ 
tor.  It  is  this  value  of  the  predicted  concentration  variance  (or,  equivalently,  mean-square 
concentration)  that  should  be  used  in  Eq.  (32)  in  the  specification  of  f{c^\S). 

It  should  be  noted  that  the  averaging  imposed  by  a  slower  response  detector  (e.g.,  with 
response  time  greater  than  a  few  seconds)  depends  critically  on  the  concentration  integral 
time  scale  Tc-  For  example,  at  short  range  in  unstable  conditions,  the  highest  concentrations 
have  very  short  durations  (and,  hence,  small  integral  time  scales),  and  consequently  the 
time  averaging  imposed  by  a  detector  may  be  significant  in  this  case  because  concentration 
fluctuations  shorter  than  about  1  s  are  averaged  in  most  chemical  agent  detectors  (e.g., 
CAM  and  AP4C  which  have  typical  averaging  times  of  5-10  s).  However,  under  stable 
conditions,  the  peak  concentrations  have  much  longer  durations  and  are  also  higher  relative 
to  the  mean  concentration.  Under  these  conditions,  the  integral  time  scales  of  concentration 
are  longer  and,  as  a  consequence,  the  concentration  fluctuations  would  not  be  smoothed  out 
by  a  chemical  agent  detector  to  the  same  extent.  The  same  situation  applies  to  dispersion 
in  a  built-up  (urban)  environment,  where  vigorous  mixing  of  plume  elements  in  canopy 
turbulence  results  in  a  dramatic  increase  in  the  integral  time  and  length  scales  here  [17]. 
More  specihcally,  the  enhanced  and  rapid  mixing  of  plume  material  in  a  built-up  area  merges 
small-scale  plume  concentration  structures,  leading  to  a  reduction  in  the  frequency  and 
number  of  low-amplitude  events  and  high  amplitude  bursts  relative  to  that  in  a  comparable 
open-terrain  plume,  the  net  effect  of  which  is  to  increase  the  concentration  fluctuation 
integral  time  scale  Tc- 

7  Conclusions 


This  report  examines  the  form  of  the  concentration  PDF  for  plumes  dispersing  in  various 
obstacle  arrays  (built-up  environment).  The  form  of  the  concentration  PDF  was  studied 
by  examining  the  relationships  between  various  normalized  higher-order  concentration  mo¬ 
ments  using  a  large  and  comprehensive  data  set  obtained  from  a  series  of  measurements 
of  plume  dispersion  in  a  number  of  obstacle  arrays  made  in  a  boundary-layer  water  chan¬ 
nel  using  high-resolution  laser-induced  fluorescence.  For  the  present  study,  we  extracted  a 
large  number  of  concentration  time  series  from  a  number  of  different  experiments  involving 
various  obstacle  arrays  (e.g.,  various  arrays  of  cubical  and  non-cubical  obstacles  in  aligned 
and  staggered  arrangements  with  uniform,  random,  and  alternating  heights). 

The  key  conclusions  revealed  by  this  analysis  can  be  summarized  as  follows.  A  remarkably 
robust  feature  of  all  the  concentration  data  was  the  observed  collapse  of  the  third-  and 
fourth-order  normalized  concentration  moments  on  the  second-order  normalized  concen- 
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tration  moment  and  of  the  observed  {S,  K)  collapse  to  a  series  of  “universal”  curves.  The 
collapse  of  the  various  concentration  moments  here  was  exactly  the  same  as  that  observed  in 
the  CONFLUX  concentration  data  for  open-terrain  plumes.  This  remarkable  collapse  sug¬ 
gests  that  the  concentration  PDF  of  plumes  dispersing  in  either  a  built-up  or  open-terrain 
environment  can  be  described  adequately  by  at  most  two  parameters  (namely,  a  location 
parameter  which  can  be  chosen  to  be  the  mean  concentration  and  a  scale  parameter  which 
can  be  chosen  to  be  the  root-mean-square  concentration  or,  equivalently,  the  concentration 
standard  deviation).  Finally,  the  observed  relationships  between  the  various  higher-order 
normalized  concentration  moments  were  found  to  be  adequately  modeled  (predicted)  using 
either  a  clipped-gamma  distribution  or  the  simpler  intermittent  exponential  distribution 
(which  is  a  special  case  of  the  clipped-gamma  distribution).  However,  the  general  shape 
of  the  observed  concentration  probability  distribution  is  well  modeled  using  the  clipped- 
gamma  distribution,  but  the  intermittent  exponential  distribution  provides  a  rather  poor 
conformance  with  the  observed  concentration.  Indeed,  the  clipped-gamma  distribution  cap¬ 
tures  correctly  the  form  of  both  the  lower  and  upper  tails  of  the  measured  concentration 
probability  distribution  for  array  plumes. 

The  scientific  case  for  using  probabilistic/statistical  methods  to  assess  the  consequences  of  a 
toxic  release,  rather  than  traditional  deterministic  methods,  is  overwhelming.  In  this  paper, 
the  implication  of  the  form  of  the  concentration  PDF,  for  the  formulation  of  a  probabilistic 
model  for  the  response  of  a  chemical  agent  detector,  is  investigated  in  the  context  of  the 
source  reconstruction  problem.  To  this  purpose,  the  likelihood  function  of  a  chemical  agent 
detector  (e.g.,  CAM,  AP4C),  which  only  reports  short  time-averaged  concentration  levels 
discretized  on  an  A^fe-level  bar  scale,  has  been  derived.  The  close  connection  between  the 
likelihood  function  and  the  concentration  PDF  has  been  investigated  in  detail. 

The  likelihood  function,  which  characterises  the  likelihood  that  a  chemical  agent  detector 
will  give  a  certain  bar  level  reading  for  a  hazardous  gas  release  from  a  given  source  distri¬ 
bution  5,  was  derived.  The  derivation  explicitly  accounts  for  (1)  model/input  uncertainties 
in  prediction  of  C  required  in  the  definition  of  the  concentration  PDF;  (2)  model/input 
uncertainties  in  prediction  of  required  in  the  definition  of  the  concentration  PDF;  (3) 
measurement  noise  in  the  detector;  and,  (4)  time-averaging  imposed  by  the  finite  response 
of  the  detector.  This  likelihood  function  can  be  used  directly  in  Bayesian  source  inversion 
methodologies  [27]  that  have  been  developed  to  determine  unknown  source  characteristics 
(e.g.,  location,  source  strength,  time  of  release,  etc.)  when  given  a  finite  number  of  noisy 
concentration  measurements  made  by  an  array  of  detectors. 
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